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Abstract 


Quantitative  processing  of  signal  space  data  to  pro- 
vide transmission  line  analysis  was  developed  in  support  of 
work  being  done  at  the  United  States  Air  Force  Rome  Air 
Development  Center  (RADC)  laboratories.  The  RADC  engineers 
were  working  with  digital  data  modems  and  they  were  quali- 
tatively analyzing  transmission  line  perturbations  through 
observation  of  an  oscilloscope  display  of  the  signal  space 
representation.  The  display  was  being  generated  by  an 
optional  circuit  available  with  the  modem,  but  the  oscillo- 
scope display  was  not  capable  of  providing  adequate  quanti- 
tative perturbation  information.  However,  the  signal  space 
data  did  contain  all  of  the  information  required  for  a 
complete  quantitative  analysis  of  xhe  transmission  line 
perturbations.  Therefore,  a system  using  the  same  data 
potentially  could  be  devised  to  perform  the  analysis. 

A"  signal  space  data  analysis  system  design  was  devel- 
oped in  this  thesis,  which  would  be  capable  of  providing 
three  displays  of  the  signal  space  representation,  the 
amplitude  and  phase  deviations,  and  the  frequency  spectrum 
of  the  phase  deviations  or  phase  jitter.  The  system  would 
consist  of  the  digital  modem,  a minicomputer,  and  an 
interface  device  between  the  modem  and  minicomputer  which 
would  be  capable  of  interrupting  the  minicomputer.  A 
FORTRAN  program  was  also  developed  and  run  in  a simulation 
effort  which  transformed  the  xy  data  of  the  signal  space 
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into  a form  of  data  vvhich  could  be  plotted  on  the  three 
displays.  Though  the  system  was  not  actually  constructed, 
it  was  successfully  simulated  such  that  it  would  merit 
implementation. 
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I . Introduction 


Through  analysis  of  the  changes  in  the  coordinates  of  a 
received  data  communications  signal  in  signal  space,  one  can 
obtain  several  characteristics  of  the  transmission  line  or 
channel  over  which  the  signal  was  transmitted.  This  is  pos- 
sible because  channel  perturbations  such  as  noise,  phase 
jitter,  hits,  harmonic  distortion,  and  line  outages  will 
affect  the  signal  space  representation  of  a data  signal  in 
unique  and  consistant  manners.  Also,  since  most  systems  use 
data  signal  point  randomizers,  the  signal  space  pattern  is 
independent  of  the  data.  The  signal  space  patterns  of  an 
active  line  can  be  analyzed  without  knowing  the  data  content 
and  without  interfering  with  the  data  transmission.  Thus, 
the  line  perturbations  can  be  determined  without  removing 
the  line  from  service  as  must  be  done  to  perform  conven- 
tional transmission  line  measurements.  However,  the  signal 
space  analysis  must  be  performed  in  real  tine  to  completely 
obtain  the  line  characteristics,  and  the  required  speed  has 
not  been  practical  until  the  low  cost  development  of  mini- 
computers. Practical  real  time  analysis  of  the  signal  space 
for  data  rates  up  to  9600  bits  per  second  has  been  possible 
since  1976  within  the  International  Business  Machines  (IBM) 
corporation  (Ref.  l),  but  a method  of  performing  the  anal- 
ysis has  not  been  made  available  outside  the  IBM  compainy. 

This  thesis  was  developed  to  provide  a system  of- anal- 
ysis which  could  be  used  to  analyze  in  real  time  the  signal 
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space  representation  of  a data  signal  and  which  could  pro- 
vide quantitative  line  characteristic  outputs  much  like  the 
real  time  outputs  of  the  IBM  system.  The  signal  space  coor- 
dinates could  be  taken  from  a digital  modem  receiver  in 
binary  form  and  fed  to  a computer  for  comparison  with  the 
known  transmitted  signal  space  coordinates  (since  the  signal 
space  coordinates  are  independent  of  the  data,  the  data 
symbols  do  not  need  to  be  determined  to  perform  the  anal- 
ysis). The  difference  between  the  transmitted  and  received 
signal  space  coordinates  could  then  be  analyzed  to  determine 
the  characteristics  of  perturbations  on  the  transmission 
line.  The  signal  space  data  could  be  taken  from  the  modem 
via  an  external  transfer  register  which  could  tv  read  by  the 
computer  when  an  interrupt  flag  is  set  by  the  modem. 

A FORTRAN  program  which  is  capable  of  analyzing  the 
signal  space  data  and  providing  a quantitative  output  was 
developed  for  this  thesis  but  it  must  be  expanded  to  work 
optimally  on  the  particular  minicomputer  which  will  be  util- 
ized in  the  laboratory.  The  FORTRAN  program  was  developed 
for  use  on  a multiuser  computer  system  which  did  not  have  an 
interrupt  capability,  therefore  the  program  does  not  contain 
the  interrupt  handling  routines  required  to  read  the  input 
data  and  to  control  the  output  displays.  The  required 
program  expansions  are  discussed  in  Section  III. 

The  output  created  by  the  FORTRAN  program  consists  of 
real  time  printer  displays  of  signal  space,  of  the  amplitude 
and  phase  deviation  for  each  data  point,  and  of  the 
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frequency  spectrum  for  phase  deviations.  Though  the  program 
is  currently  able  to  provide  only  real  time  displays  of 
transmission  line  perturbations,  the  analysis  system  would 
be  greatly  enhanced  if  the  capabilities  of  averaging,  creat- 
ing, storing,  and  displaying  a time  condensed  history  of 
line  data  are  added  to  the  program.  The  structure  of  the 
program  is  designed  so  that  these  capabilities  can  be  added 
to  the  existing  instructions  rather  than  reaccomplishing  the 
entire  program. 

The  thesis  is  organized  with  sections  on  problem  back- 
ground and  definition,  solution,  results  and  recommenda- 
tions. In  the  section  on  background  and  definition  of  the 
requirement  is  a description  of  the  origin  of  the  require- 
ment and  the  practical  aspects  of  equipment  which  affected 
the  system  design.  In  the  next  section  is  a description  of 
the  development  of  the  system  structure  and  program  design, 
followed  by  a separate  section  which  centers  on  discussion 
of  program  results.  The  final  section  contains  recommenda- 
tions for  further  development  and  implementation  of  the 
system. 


II.  Background  and  Definition  of  Requirement 

The  need  for  quantivative  processing  of  signal  space 
data  to  provide  transmission  line  analysis  arose  from  work 
being  done  at  the  United  States  Air  Force  Rome  Air  Develop- 
ment Center's  (RADC)  telecommunications  digital  laboratory. 
The  telecommunications  engineers  were  doing  development 
work  in  the  laboratory  with  digital  modems  which  were  pro- 
viding a qualitative  oscilloscope  picture  of  the  real  time, 
received  signal  space  representation.  Ideally,  the  signal 
space  picture  could  show  the  condition  of  the  transmission 
line,  and  it  could  be  used  to  detect  several  transmission 
line  problems.  However,  the  qualitative  information  pro- 
vided by  the  oscilloscope  could  only  be  utilized  to  give  a 
rough  indication  of  transmission  line  quality  since  it  only 
provided  an  average  of  data  which  the  human  eye  could  detect. 
Even  a trained  person  watching  the  scope  could  not  tell 
whether  or  not  the  line  was  about  to  experience  a lineout 
(loss  of  transmission)  due  to  a gradual  degradation  in  line 
quality.  Thus  it  was  decided  that  a quantitative  method  of 
analyzing  the  signal  space  pattern  should  be  developed. 

Modem  Parameters  (Ref.  2) 

The  modems  being  used  in  the  RADC  laboratory  are  prima- 
rily Codex  LSI  9600  modems  and  are  on  consignment  to  RADC. 
They  are  full  duplex  digital  modems  with  quadrature  phase 
shift  and  amplitude  modulation  at  a center  frequency  of 
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Figure  1.  Signal  Space  Patterns  (Eye  Patterns)  for  Different 
Data  Modes  on  a Codex  96OO  LSI  Modem  (Ref.  2) 


1706  Hz  and  with  a bandwidth  of  2400  Hz.  They  are  designed 
to  work  on  C2  Conditioned  4 Wire  telephone  circuits.  The 
modems  are  capable  of  operating  at  data  rates  of  96OO  BPS 
(bits  per  second)  at  2400  baud,  or  7200  BPS  at  2400  baud,  or 
4800  BPS  at  I600  baud,  or  4800  BPS  at  2400  baud  all  of  which 
are  switch  selectable  on  the  modem  control  panel.  The  cor- 
responding signal  space  pattern  for  these  different  data 
rates  contain  16  points  for  the  highest  data  rate  of  96OO 
BPS,  8 points  for  the  next  two  rates  and  4 points  for  the 
lowest  rate  (Fig.  l).  The  modems  contain  equalizers  for 
adjusting  the  delay  and  for  equalizing  the  amplitude  between 
the  two  quadrature  components  during  each  baud  time.  They 
also  contain  data  randomizers  which  rotate  the  different 
signal  space  points  for  each  baud  so  that  the  data  is 
uniformly  distributed  among  all  the  signal  space  points. 
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The  signal  space  pattern  is  generated  by  a Codex  system  as 
an  optional  circuit  inside  the  modems  or  as  an  optional 
external  plug-in  circuit.  The  digital  signal  which  feeds 
the  generator  is  taken  from  inside  the  modem  as  x and  y 
values  of  signal  space  after  demodulation  and  delay  and 
quadrature  equalization  (Fig.  2). 

Operation  of  Modem  (Ref.  2).  The  analog  received  data 
signal  is  converted  to  a serial  digital  data  string  in  the 
sample  and  hold  circuit  prior  to  the  demodulation  section 
(digital  multiply)  (Fig.  2).  Then  the  serial  data  is  de- 
modulated using  digital  filter  algorithms.  The  adaptive 
equalizer  adjusts  each  point  of  signal  space  information 
(equivalent  to  each  baud)  for  delay  and  amplitude  equali- 
zation between  the  two  quadrature  components  as  it  is  fed 
into  the  xy  accumulator  (the  x and  y values  for  the  signal 
space  point  are  multiplexed  into  a single  word  of  informa- 
tion) . Word  synchronization  is  supplied  by  another  circuit 
within  the  modem  and  signals  the  accumulator  to  feed  the 
signal  space  word  into  the  data  decision  logic  at  the  appro- 
priate time.  The  data  decision  logic  simply  demultiplexes 
the  xy  values,  determines  the  closest  preset  signal  space 
point  to  the  received  signal  space  point,  and  assigns  the 
appropriate  data  value  as  determined  by  the  randomizer 
circuit  (not  shown  in  Fig.  2).  The  data  decision  logic  is 
'simply  a minimum  distance  decision  system. 

The  signal  space  point  coordinates  are  stored  in  ROM 
inside  the  modem  and  are  set  at  the  factory.  As  different 
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data  rate  modes  are  selected,  the  signal  space  coordinates 
do  not  change,  only  the  number  of  points  utilized  for  trans- 
mission are  changed,  and  only  the  randomizer  makes  the 
appropriate  changes.  As  shown  in  Figure  1 certain  points 
are  eliminated  as  the  data  rate  is  slowed  down.  Also,  the 
number  of  bits  represented  by  each  signal  point  is  changed 
from  4 to  3 to  2 depending  on  the  data  mode  selected. 
However,  this  again  only  affects  the  randomizer  mechanism 
and  not  the  signal  space  mechanism  since  the  same  signal 
space  coordinates  are  used  in  each  mode.  For  transmission, 
the  xy  coordinates  of  a point  are  put  into  two  words  of  five 
bits  per  word.  Then  the  randomizer  selects  a different 
signal  point  for  each  baud  of  information  and  the  modem 
transforms  and  transmits  the  two  words  of  signal  space  coor- 
dinates as  quadrature  sin  and  cos  frequencies  with  magnitude 
and  phase  information.  At  the  receiver,  two  words  of  eight 
bits  each  are  used  to  describe  the  signal  space  xy  coor- 
dinates' of  the  received  signal.  These  coordinates  are  then 
compared  with  the  coordinates  stored  in  ROM  to  determine  the 
closest  transmitted  signal  space  point.  The  two  words  of 
eight  bits  each  are  also  used  to  generate  the  signals  for 
the  oscilloscope  signal  space  picture. 

Operation  of  Eve  Pattern  Generator  (Ref.  2).  The  two 
Codex  options  for  generating  the  oscilloscope  signal  space 
picture  both  use  the  same  circuit  configuration.  The  in- 
ternal option  is  merely  included  on  a plug-in  card  inside 
each  modem  case  and  the  external  option  puts  the  circuit  in 
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a small  case  that  can  be  plugged  into  the  external  jacks  on 
any  Codex  LSI  96OO  modem.  The  Codex  company  calls  the 
options  "Eye  Pattern  Generators".  The  generators  consist  of 
a demultiplexer  circuit,  two  serial  to  parallel  converters, 
two  digital  to  analog  converters,  and  buffers  and  amplifiers 
for  interface  with  the  modem  and  oscilloscope  (Fig.  3).  The 
xy  serial  signal  space  data  is  fed  into  the  Eye  Pattern 
Generators  from  the  modem  xy  accumulator  (Fig.  2)  and  the 
word  synchronization  for  the  generator  is  supplied  from  the 
modem  word  synchronization  signal.  The  demultiplexer  breaks 
the  serial  xy  data  word  into  two  words  for  x and  y of  eight 
bits  each  and  strings  them  into  the  serial  to  parallel 
registers.  Then,  the  word  synch  signal  from  the  demulti- 
plexer signals  the  converters  to  simultaneously  dump  the  8 
bits  of  X and  8 bits  of  y coordinates  into  the  digital  to 
analog  converters.  These  converters  in  turn  generate  a 
voltage  level  which  can  be  amplified  to  drive  an  ordinary 
oscilloscope  to  generate  the  signal  space  pattern  or  eye 
pattern. 

Eve  Pattern  Characteristics  (Ref.  2) 

Since  the  signal  space  points  are  independent  of  the 
data,  the  line  perturbations  will  be  contained  in  the  signal 
fed  to  the  oscilloscope  and  will  not  be  affected  by  nor  in- 
'terfere  with  the  data  being  transmitted  over  the  line.  The 
only  perturbations  which  will  not  show  up  on  the  oscillo- 
scope are  differences  in  quadrature  amplitude  distortions 
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and  delay  because  the  adaptive  equalizer  corrects  them  prior 
to  the  point  where  the  data  is  fed  to  the  Eye  Pattern  Gener- 
ator. Pure  line  perturbations  of  lineout,  normal  gaussian 
noise,  impulse  noise,  phase  jitter,  phase  hits,  and  ampli- 
tude hits  will  show  up  on  the  oscilloscope  signal  space 
pattern  as  depicted  in  Figure  4 for  7200  BPS  at  2400  baud  or 
4800  BPS  at  1600  baud. 

Lineout . Lineout  occurs  when  the  signal  loses  synchro- 
nization and  the  transmitter  ceases  to  transmit.  Therefore, 
the  only  received  signal  is  the  noise  on  the  line  which 
contains  various  phases  and  small  amplitudes  which  center 
around  the  origin  of  the  signal  space. 

Normal  Gaussian  Noise . Normal  gaussian  noise  adds  and 
subtracts  from  the  transmitted  signal  amplitude  and  phase. 
Therefore,  the  received  signal  points  will  have  a random 
scattering  about  their  respective  transmitted  signal  space 
points . 

Impulse  Noise ■ Impulse  noise  is  very  narrow  in  time 
and  consequently  only  affects  one  or  two  received  signal 
space  points.  Thus,  the  impulse  noise  appears  as  a momen- 
tary single  spot  on  the  oscilloscope  far  away  from  the 
general  cluster  of  received  points. 

Phase  Jitter.  Phase  jitter  is  the  oscillation  of  the 
sequential  transmitted  phase  points  about  their  respective 
reference  phase  points.  The  oscillation  v/ill  simply  cause  a 
rotational  rocking  motion  of  the  entire  signal  space  pattern 
as  it  is  observed  on  the  oscilloscope. 
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Figure  4.  Oscilloscope  Digital  Eye  Pattern  (Signal  Space 
Pattern)  with  Pure  Signal  and  with  Several  Pure 
Line  Perturbations  (Ref.  2) 
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Phase  Hit . A phase  hit  is  caused  by  a constant  phase 
angle  being  added  to  or  subtracted  from  each  transmitted 
point  for  a short  period  of  time.  The  affect  on  the  oscil- 
loscope is  that  the  signal  space  pattern  rotates  by  the 
amount  of  the  hit,  then  rotates  back  to  normal. 

Amplitude  Hit . An  amplitude  hit  is  the  same  as  a phase 
hit  except  that  the  constant  deviation  affects  the  ampli- 
tude. The  affect  on  the  oscilloscope  is  that  the  signal 
space  pattern  seems  to  expand  or  shrink  about  the  signal 
space  origin  then  return  to  normal. 

Limitations  of  the  Eve  Pattern.  Though  the  pure  cases 
of  perturbations  are  easy  to  see  and  detect  on  the  oscillo- 
scope, they  cannot  be  quantitatively  analyzed  by  an  ob- 
server. Also,  when  several  different  perturbations  occur 
simultaneously,  the  types  of  perturbations  cannot  be  deter- 
mined, and  when  the  perturbations  happen  at  a high  rate,  the 
human  eye  cannot  follow  the  changes.  For  instance,  if  the 
line  is  noisy  and  a small  phase  jitter  is  present,  the  phase 
jitter  will  not  be  seen  on  the  oscilloscope  because  the 
phase  jitter  will  be  covered  up  by  the  noise.  If  impulse 
noise  strikes  only  very  sporadically^  it  will  not  even  be 
seen  on  the  scope  since  one  or  two  points  will  appear  for 
only  0.4  milliseconds  at  2400  baud.  Therefore,  the  oscillo- 
scope signal  space  picture  can  only  be  used  to  roughly 
detect  line  perturbations  and  it  only  indicates  that  a 
general  problem  is  occuring  on  the  transmission  line. 
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statement  of  Requirement 

The  primary  problem  is  that  the  oscilloscope  signal 
space  picture  does  not  provide  enough  quantitative  infor- 
mation to  fully  analyze  the  transmission  line  condition. 
However,  the  information  provided  to  the  oscilloscope  does 
contain  most  of  the  line  perturbation  information  which  is 
required  for  a complete  analysis  of  the  line.  Also,  since 
the  modems  are  only  on  consignment  to  the  RADC  laboratory, 
no  internal  modifications  may  be  made  to  the  modems  and 
therefore  only  the  information  available  at  the  output 
points  may  be  utilized.  Therefore,  the  problem  is  to  devise 
a method  of  quantitative  analysis  which  uses  either  the 
digital  eye  pattern  signal  or  else  the  analog  eye  pattern 
signal  to  provide  quantitative  information  about  trans- 
^ mission  line  perturbations. 
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III.  Development  of  Solution 

The  development  of  a system  which  could  quantitatively 
analyze  the  signal  space  data  to  obtain  a transmission  line 
analysis  was  straight  forward.  Interviews  with  the  RADC , 
Codex  and  IBM  engineers,  and  research  of  literature  led 
directly  to  a formulation  of  the  system  design  and  to  the 
decision  of  analyzing  the  digital  signal  space  data  from  the 
modem  with  a digital  computer.  Since  the  basic  concept  of 
analyzing  signal  space  data  for  transmission  line  analysis 
could  apply  to  any  digital  modem  using  minimum  distance 
decision  logic,  it  was  also  decided  to  keep  the  program 
flexible  so  that  it  could  be  adapted  to  different  modems 
and  computers.  Therefore,  the  FORTRAN  program  was  developed 
' such  that  it  contained  functionally  structured  elements  of 
routines  that  could  be  adapted  or  expanded  without  major 
modification  to  the  entire  program. 

Formulation  of  the  System  Design 

During  the  initial  interviews  with  the  RADC  engineers 
and  Codex  engineers,  it  was  learned  that  the  IBM  company  had 
published  an  article  in  late  19(6  (Ref.  l)  which  described  a 
transmission  line  monitor  concept  utilizing  signal  space 
data.  An  interview  with  the  IBM  engineers  was  arranged  and 
their  LQM  (Transmission  Line  Quality  Monitor)  system  was 
demonstrated.  Their  LQM  system  demonstration  was  impressive 
in  that  it  provided  a great  deal  of  real  time  analysis 


15 


outputs  and  of  condensed  historical  numerical  outputs  on  the 
transmission  line  quality  and  performance.  However  the  LQM 
system  was  extremely  complex  and  it  was  strictly  an  internal 
IBM  system,  thus  details  of  the  system  operation  were  not 
available  outside  the  IBM  company.  But,  since  a workable 
system  of  digital  analysis  of  signal  space  data  for  line 
analysis  had  not  been  demonstrated,  it  was  decided  that  a 
simpler  system  using  the  same  concepts  to  provide  only  real 
time  analysis  outputs  was  within  the  scope  of  a thesis  and 
that  such  a syatem  could  be  developed  for  the  RADC  labora- 
tory requirements. 

System  Concepts . The  basic  concept  of  the  system  de- 
sign is  that  since  most  of  the  transmission  line  character- 
istics are  contained  in  the  two  x and  y words  which  describe 
- each  received  data  point  in  the  signal  space,  the  analysis 
system  need  only  be  comprised  of  a device  for  extracting  the 
digital  signal  space  data  from  the  modem  and  a device  for 
analyzing  and  displaying  the  processed  data  to  show  the  line 
characteristics.  This  concept  meets  the  constraints  of  not 
modifying  the  modems  and  of  using  raw  data  available  at  the 
modem.  , 

System  Hardware . The  data  extraction  device  would  be 
similar  to  the  Codex  optional  external  Eye  Pattern  Genera- 
tor, and  the  analysis  and  display  device  could  consist  of 
any  minicomputer  comparable  with  the  Digital  Equipment 
Corporation's  PDP-ll/45  minicomputer  which  is  available  in 
the  RADC  laboratories  (the  IBM  system  operates  on  the  small 
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IBM  System-7)  (Ref.  l). 

The  data  extraction  mechanism  would  have  to  be  able  to 
input  the  serial  digital  signal  space  data  from  the  modem 
output  port  for  the  external  Eye  Pattern  Generator,  convert 
it  to  two  words  of  signal  space  coordinates,  and  feed  the 
coordinates  to  the  minicomputer  for  processing.  As  in  the 
external  Eye  Pattern  Generator  (Fig.  3) . the  data  would  have 
to  be  demultiplexed  and  converted  to  two  words  of  parallel 
data.  As  shown  in  Figure  5 "the  extraction  mechanism  would 
contain  input  buffers,  a demultiplexer,  two  serial  to 
parallel  converters,  plus  output  buffers  and  handshaking 
circuits  to  match  with  the  computer  interrupt  protocol.  A 
det'lled  design  was  not  developed  in  this  thesis  because  the 
working  design  would  incorporate  several  mechanisms  already 
. available  in  the  RADC  laboratories  which  could  be  more 

efficiently  interfaced  by  the  lab  technicians  who  are  famil- 
iar with  them. 

The  minicomputer  which  would  be  utilized  for  signal 
space  analysis  and  display  would  have  to  have  an  interrupt 
capability  or  possibly  a direct  memory  access  capability 
(DMA).  The  interrupt  or  DMA  capability  would  allow  the 
signal  space  coordinates  to  be  fed  into  the  computer  as  they 
would  be  received  in  the  data  extraction/interface  device 
one  point  at  a time.  At  the  highest  data  rate  of  2400  baud, 
there  would  be  4l6  microseconds  between  data  points.  There- 
fore the  interrupt  handler  of  the  minicomputer  would  have  to 
complete  the  servicing  routine  in  that  time  frame  or 
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subsequent  data  points  would  be  lost.  The  worst  case  time 
for  a PDP-ll/45  minicomputer  to  acknowledge  and  branch  to 
the  first  instruction  of  an  interrupt  servicing  routine 
would  be  1?.19  microseconds.  The  basic  move  of  the  input 
data  from  the  external  transfer  register  to  a location  in 
the  memory  would  only  be  required  once  and  the  worst  case 
time  would  be  4.28  microseconds.  Therefore,  there  would  be 
394.53  microseconds  left  over  for  the  basic  operations  of 
the  interrupt  handler  routine  and  indexing  of  the  input 
buffer  stack  pointer  (Ref.  5)*  These  time  frames  would  be 
adequate  for  the  simpler  current  data  analysis  system  devel- 
oped in  this  thesis  since  the  subsequant  data  points  after 
the  initial  filling  of  the  input  buffer  would  not  be  needed. 
However,  for  an  expanded  analysis  system,  these  time  frames 
would  be  critical  and  require  that  more  elaborate  input 
buffering  techniques  be  developed.  DMA  would  be  consider- 
ably faster,  but  a more  complex  interface  would  be  required 
between  the  modem  and  computer  (Ref.  6). 

The  output  of  the  analysis  system  would  consist  of 
graphical  displays  showing  the  desired  quantitative  informa- 
tion on  the  line  perturbations.  These  displays  could  be 
shown  on  a video  screen  or  printed  out  in  hard  copy.  The 
FORTRAN  program  included  in  this  thesis  was  designed  to 
print  the  displays  on  the  line  printer. 

Development  of  Program  Structure 

The  development  of  the  concepts  and  equipment 
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requirements  for  real  time  signal  space  analysis  were 
straight  forward,  however  the  development  of  the  computer 
program  which  would  perform  the  analysis  was  a much  more 
difficult  procedure.  The  program  had  to  be  capable  of 
converting  two  integer  values  of  x and  y into  parameters 
which  would  describe  the  transmission  line  condition  in 
standard  terms  and  it  had  to  be  done  in  such  a way  that  the 
program  could  be  expanded  at  a later  date  to  include  the 
capability  of  providing  historical  data  on  the  transmission 
line  performance.  Therefore,  it  was  decided  to  use  design 
techniques  which  were  developed  in  a previous  AFIT  thesis 
(Ref.  3)  to  define  and  develop  the  program  requirements  and 
parameters.  Also,  these  techniques  would  allow  the  program 
to  be  expanded  and  adapted  to  different  computer  systems  at 
a later  time  much  easier  than  if  just  a single  real  time 
analysis  program  were  developed  for  a particular  machine. 

The  development  of  the  program  structure  first  consisted  of 
defining  all  of  the  requirements  and  parameters  of  the  IBM- 
LQM  system  so  that  the  entire  scope  of  an  expanded  system 
could  be  included  in  the  basic  design.  Then  a structured 
programming  technique  was  applied  which  allowed  the  develop- 
ment of  the  diagrams  in  Appendix  A.  These  diagrams  were 
then  converted  into  a basic  structure  chart  as  described  in 
Reference  3*  The  basic  structure  chart  v/as  then  refined  so 
that  it  would  make  possible  the  straight  forward  programming 
of  the  real  time  analysis  or  current  data  processing  sec- 
tions of  the  overall  system.  Since  only  the  definition 
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phase  and  the  final  program  structure  were  peculiar  to  this 
thesis,  the  refinement  phases  were  not  developed  in  the 
thesis  text. 

Definition  of  Program  Requirements  and  Parameters . The 
program  input,  output,  and  control  requirements  and  param- 
eters were  developed  through  an  iterative  process  of  pro- 
gressively increasing  the  detail  of  functional  diagrams. 
Three  levels  of  these  diagrams  are  shown  in  Appendix  A,  but 
only  the  last  level  is  explained  here  as  it  contains  all  of 
the  program  requirements  and  parameters  of  a complete  sys- 
tem. 

The  data  input  would  consist  of  the  two  values  of  x and 
y which  describe  a single  point  in  the  signal  space  (for  the 
Codex  9600  LSI  modem,  these  two  words  have  eight  bits  each 
and  therefore  have  integer  values  between  plus  or  minus 
128).  The  data  input  for  a complete  system  would  be  started 
at  time  zero  and  run  continuously  with  x and  y values  being 
received  simultaneously  at  regular  rates  of  2400  or  16OO 
baud,  depending  on  the  data  rate  of  the  modem.  However,  the 
data  rate  would  not  affect  the  operation  of  the  program 
since  it  is  designed  for  use  at  the  highest  rate.  It  should 
be  noted  here,  that  the  current  data  processing  elements 
would  not  need  to  operate  continuously  if  they  were  all  that 
were  to  be  developed  for  the  analysis  system.  The  data 
input  would  simply  be  started  when  the  operator  requested, 
ajid  the  input  routine  would  only  run  until  an  input  buffer 
was  filled. 
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The  output  of  the  program  would  consist  of  displays 
which  show  standard  parameters  of  transmission  lines.  It 
was  decided  that  the  displays  developed  by  the  IBM-LQM 
system  provided  all  of  the  needed  transmission  line  param- 
eters, so  the  output  displays  of  this  program  structure  were 
patterened  after  the  IBM  displays.  The  displays  would  con- 
sist of  real  time  displays  showing  the  signal  space,  the 
amplitude  and  phase  deviations,  and  the  phase  deviation 
frequency  spectrum.  There  would  also  be  displays  which 
would  show  summaries  of  five  minutes,  one  hour,  twenty  four 
hours  of  a single  transmission  line  (the  IBM-LQM  system  can 
monitor  up  to  seven  lines),  and  a complete  summary  of  all 
lines  for  a twenty  four  hour  period.  These  summaries  would 
contain  num^ical  averages  and  standard  deviations  for 
amplitude  deviation,  phase  deviation,  phase  jitter,  random 
noise,  phase  dispersion,  and  line  quality.  The  summaries 
would  also  contain  the  number  of  severe  or  medium  phase  and 
amplitude  hits  for  the  particular  period  of  the  summary. 
Though  the  summaries  would  add  a great  deal  of  information 
capability  to  the  analysis  system,  they  would  also  require 
several  thousand  man  hours  to  be  developed  (Ref.  1).  There- 
fore they  are  only  described  here  to  show  that  they  were 
considered  in  the  basic  design,  but  not  with  the  intent  of 
being  developed  in  this  thesis.  The  displays  could  either 
be  shown  on  a video  screen  or  be  printed  on  a hard  copy. 

To  convert  the  input  from  xy  data  to  quantitative 
output  displays,  several  operations  would  have  to  be 
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performed  on  each  data  point  and  on  different  groupings  of 
data  points  as  they  would  be  sequentially  received.  For  the 
current  data  displays,  the  xy  coordinates  would  need  to  be 
converted  to  polar  coordinates,  and  the  deviations  from  the 
signal  space  reference  points  or  target  points  would  need  to 
be  determined  in  polar  coordinates.  Then  the  phase  devia- 
tions would  have  to  be  fast  Fourier  transformed  to  determine 
the  phase  deviation  frequency  spectrum  or  phase  jitter 
frequencies.  Finally,  the  data  would  need  to  be  arrayed  so 
that  it  could  be  displayed  in  graphibal  form.  If  the 
summary  displays  and  historical  capabilities  were  to  be 
added  to  the  system,  then  the  current  processed  data,  the 
original  raw  xy  data,  plus  additional  statistical  processed 
data  would  have  to  be  processed,  arrayed,  and  stored.  To 
obtain  the  final  configuration  of  the  functional  diagrams  as 
shown  in  Appendix  A and  of  the  program  structure  diagrams  as 
shown  in  Figure  6 through  10,  several  iterations  of  design 
were  required.  These  iterations  included  the  use  of  bubble 
charts  and  several  different  program  structure  charts,  how- 
ever, it  was  decided  that  complete  explanations  of  the 
iterations  were  not  required  to  understand  the  final  program 
structure . 

Program  Structure  Diagrams . The  top  level  program 
structure  diagram  (Fig.  6)  is  arranged  with  the  input  or 
efferent  element  on  the  left  and  with  the  output  or  afferent 
element  on  the  right  so  that  data  flow  and  order  of  process- 
ing operations  would  tend  to  be  left  to  right.  Each  of  the 


23 


r«trforrf\ 

Trar\&n\iiSlori 

Line  /\r\a\ysis 


TCsrath’ 


mm 


/ 


i.sf 

kz8  V/o 


In  ifiaiiitc  ftnd  \Ff£>ccs>s  Process  \Crsa1& 

Start  Arai/si5  Cwrrer^  Data  Rzcafd  DaT&  Oiipts?tys 

<uh.STRTf\N  Suk.ClARODAT  Sub^.RECROAT  Suh'.  DISPLAY 


P/Tererlt 


Tranifor  mention 


A'fPe  ra.r\i 


Figure  6.  Top  Level  of  Program  Structure  Diagram 


# DESCRIPTION 


FROM 


TO 


1 . 

Initialization  Parameters 

ENTPARM 

STRTAN 

STRTAN 

TXLNAN 

TXLNAN 

CURRDAT 

RECRDAT 

2. 

Constants  for  Target 

ENTPARM 

STRTAN 

Coordinates  and  Timing 

STRTAN 

TXLNAN 

Control 

TXLNAN 

CURRDAT 

RECRDAT 

DISPLAY 

CURRDAT 

APDEV 

RECRDAT 

COMPRES 

HISTORY 

DISPLAY 

SELDAT 

S ELD AT 

SLCURR 

SLFMSM 

SLLNSM 

SLHRSM 

SLDASM 

3. 

External  Control  Data 

CONINT 

WAIT 

(Such  as  Display  Selection 

WAIT 

STRTAN 

Requests ) 

STRTAN 

TXLNAN 

TXLNAN 

DISPLAY 

DISPLAY-  — 

SELDAT 

SELDIS 

• 

SELDAT 

SLCURR 

SLFMSM 

SLLNSM 

SLHRSM 

SLDASM 

SELDIS 

DISXY 

DISAPD 

DISPDFS 

DISFMSM 

DISLNSM 

DISHRSM 

DISDASM 

Table  1.  Data  and  Control  Identifiers  for  Program  Structure 
Diagrams  (Figure  6 through  Figure  10) 
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# DESCRIPTION  FROM  TO 


4.  XY  Coordinates  of  Single  DATINT WAIT 

Received  Data  Point  WAIT DATRDY 


5.  Full  Buffer  of  120  Points  DATRDY STRTAN 

of  XY  Data  Coordinates  STRTAN TXLNAN 

TXLNAN CURRDAT 

RECRDAT 

DISPLAY 

CURRDAT XYTOPOL 

RECRDAT COMPRES 

DISPLAY SELDAT 

SELDAT 3LCURR 


6.  Amplitude  and  Phase  Data  XYTOPOL CURRDAT 

Coordinates  for  120  CURRDAT APDEV 

Received  Data  Points  TXLNAN 

TXLNAN RECRDAT 

DISPLAY 

RECRDAT COMPRES 

DISPLAY SELDAT 

SELDAT SLCURR 


7.  Amplitude  and  Phase  APDEV CURRDAT 

Deviation  Data  of  120  CURRDAT FRQSPEC 

Received  Points  from  the  TXLNAN 

Nearest  Target  Point  TXLNAN RECRDAT 

DISPLAY 

RECRDAT COMPRES 

DISPLAY SELDAT 

SELDAT SLCURR 


8.  Phase  Deviation  Frequency  FRQSPEC CURRDAT 

Data  for  120  Received  Points  CURRDAT TXLNAN 

TXLNAN RECRDAT 

DISPLAY 

RECRDAT COMPRES 

DISPLAY SELDAT 

SELDAT SLCURR 


Table  1 Continued.  Data  and  Control  Identifiers  for  Program 

Structure  Diagrams  (Figure  6 through 
Figure  10) 
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# DESCRIPTION 


FROM 


TO 


9. 

Statistical  Data 

COMPRES 

RECRDAT 

(Compressed  Data) 

RECRDAT 

HISTORY 

10. 

Catalogued  and  Line,  Date, 

HISTORY 

RECRDAT 

Time  Referenced  Summary 

RECRDAT--- 

TXLNAN 

Data 

TXLNAN 

DISPLAY 

DISPLAY— 

SELDAT 

SELDAT 

SLFMSM 

SLLNSM 

SLHRSM 

SLDASM 

11 . 

Plot  Data  for  120  Values 

SLCURR  --■ 

SELDAT 

of  Current  Data  Display 

SELDAT 

DISPLAY 

DISPLAY 

SELDIS 

SELDIS 

DISXY 

DISAPD 

DISPDFS 

12. 

Plot  Data  for  Summary 

SLFMSM 

Displays 

SLLNSM 
SLHRSM 
SLDASM 

SELDAT 

SELDAT 

DISPLAY 

DISPLAY 

SELDIS 

1 

SELDIS 

DISFMSM 

DISLNSM 

DISHRSM 

DISDASM 

Table  1 Continued.  Data  and  Control  Identifiers  for  Program 

Structiire  Diagrams  (Figure  6 through 
Figure  10) 
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levels  from  top  to  bottom  (Fig.  6 to  Fig.  10)  correspond 
roughly  with  the  levels  of  functional  diagrams  in  Appendix 
A.  This  arrangement  would  allow  the  different  functional 
operations  to  be  separated  and  ultimately  written  into 
separate  subroutines. 

Master  Control  (Fig.  6).  The  top  block  is  the  master 
control  and  it  would  determine  which  type  of  action  at  the 
next  level  would  be  performed.  The  next  level  of  four 
functions  contains  the  major  categories  of  system  activities 
which  are  data  acquisition  and  program  initialization,  real 
time  or  current  data  processing,  record  data  processing,  and 
display  output  formulation. 

Initialization  and  Analysis  Start  (Fig.  7).  The  data 
acquisition  and  program  initialization  block  would  be  the 
starting  point  of  the  analysis.  It  contains  a block  for 
initialization  of  program  constants  such  as  tne  coordinates 
of  the  target  points  or  reference  transmitted  signal  space 
points,  a block  for  an  idle  mode  when  all  other  processing 
would  have  been  completed,  and  a block  which  would  transfer 
a full  input  buffer  of  data  into  a processing  location 
within  the  computer  to  isolate  input  and  processing  actions. 
The  initialization  or  "enter  parameters"  block  v/ould  be  a 
one  time  operation  in  that  it  would  enter  all  the  needed 
constants  of  the  program  and  set  all  of  the  decision  flags 
at  the  program  start.  The  idle  or  wait  block  would  idle 
the  computer  until  an  interrupt  signal  v^as  received.  This 
block  contains  two  subdivisions  for  control  interrupts  and 
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Figure  7 • Efferent  Elements  of  Program  Structure  Diagram 
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data  interrupts.  The  control  interrupts  were  included  so 
that  the  different  displays  could  be  requested  externally. 
The  data  interrupt  would  allow  the  two  words  of  x and  y 
coordinates  to  be  read  in  from  the  external  interface  device 
and  put  into  an  input  buffer.  The  input  buffer  would  be 
necessary  in  the  expanded  system  so  that  other  processing 
actions  could  be  accomplished  while  the  buffer  was  being 
filled  by  the  data  interrupt  and  so  that  no  data  points 
would  be  lost  during  the  process.  The  last  block  in  the 
start-up  of  the  analysis  is  the  final  acquisition  of  data  or 
data  ready  block.  This  block  would  transfer  the  whole  input 
buffer  of  xy  data  into  a working  location  of  memory  so  that 
the  other  processing  operations  could  be  performed  without 
interfering  with  the  data  input. 

Current  Data  Processing  (Fig.  8).  The  current  data 
processing  block  would  prepare  the  xy  data  so  that  the 
record  data  block  could  average  and  file  the  data  and  so 
that  the  display  block  could  create  the  appropriate  displays 
of  the  current  data.  The  current  data  block  contains  three 
elements  working  in  series  which  would  convert  the  cartesian 
xy  coordinates  of  the  signal  space  to  polar  coordinates  of 
amplitude  and  phase,  calculate  the  amplitude  and  phase 
deviations  from  the  target  or  reference  transmitted  points, 
and  calculate  the  phase  deviation  frequency  spectrum. 

Record  Data  Processing  (Fig.  9).  The  record  data  block 
would  transform  the  current  data  values  of  xy  coordinates, 
amplitude  and  phase  deviations,  and  phase  deviation 
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Figure  8.  Current  Data  Elements  of  Program  Structure 
Diagram 


Figure  9.  Historical  Data  Elements  of  Program  Structure 
Diagram 
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frequencies  into  time  averages  of  deviation,  jitter,  noise, 
quality,  and  hits,  then  it  would  put  them  into  a permanent 
file.  This  block  has  two  functional  blocks  of  data  com- 
pression and  data  history  which  would  respectively  perform 
the  averaging  calculations  and  perform  the  permanent 
cataloging  or  storing  functions.  Conceivably,  these  two 
blocks  could  be  broken  down  to  several  more  levels  of  func- 
tions to  obtain  a more  readily  programmable  system.  However, 
these  blocks  would  be  highly  complex  and  would  require  an 
extensive  study  of  operations  necessary  to  calculate  the 
summary  data  values  for  the  historical  transmission  line 
analysis.  Also,  these  blocks  are  not  needed  for  the  current 
data  analysis,  so  they  are  included  only  as  expansion  ports. 

Disolav  (Fig.  10) . The  display  block  contains  two 
functional  blocks  of  which  one  would  select  the  appropriate 
data  for  a particular  display,  and  of  v/hich  the  other  would 
select  the  appropriate  display  routine  for  the  selected 
data.  The  data  selection  block  has  five  modes  of  selection 
depending  on  the  display  to  be  created.  If  the  desired 
display  was  one  of  the  three  current  data  displays,  the 
selection  block  would  gather  all  of  the  current  data  needed 
to  create  all  three  current  data  displays  since  all  three 
would  often  be  requested  for  comparison  purposes.  The  other 
selection  functions  would  be  tailored  to  each  display.  The 
display  selection  block  would  simply  call  the  appropriate 
display  subroutine  from  the  selection  of  the  seven  different 
subroutines  in  the  expanded  system,  however  only  three 
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Figure  10.  Afferent  Elements  of  Program  Structure  Diagram 
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display  subroutines  are  needed  for  a current  data  analysis 
system. 

Development  of  FORTRAN  Program 

With  use  of  the  program  structure  diagram,  the  develop- 
ment of  a functional  FORTRAN  program  was  straight  forward. 
Since  the  program  structure  contained  all  of  the  necessary 
aspects  and  parameters  of  the  system,  a complete  FORTRAN 
program  was  written  to  correspond  directly  with  each  block 
on  the  program  structure  diagram.  Though  the  resulting 
FORTRAN  program  is  not  the  fastest  for  processing  time,  nor 
most  efficient  for  memory  requirements,  it  does  contain  all 
of  the  needed  operations  required  in  the  current  data 
analysis  system  and  therefore  can  be  optimized  to  fit  a 
particular  minicomputer. 

The  first  step  in  the  FORTRAN  program  development  was 
to  assign  a subroutine  to  each  functional  block.  Then  the 
development  became  a process  of  writing  calling  routines  for 
the  top  two  levels  of  the  structure  diagram  and  of  writing 
short  functional  programs  for  each  specific  task  to  be  per- 
formed by  the  remaining  blocks.  The  top  most  block  became 
the  master  control  program  which  can  call  any  of  the  four 
blocks  or  subroutines  in  the  next  level.  These  four  sub- 
routines in  turn  can  call  the  appropriate  subroutines  under 
themselves  to  perform  the  actual  computations  or  transfers 
of  data  within  the  system.  The  details  of  the  program 
names,  decision  logic,  flag  names,  and  file  names  are 
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available  in  Appendix  B and  Appendix  C where  there  are  flow 
diagrams  and  a copy  of  the  FORTRAN  program.  Reference  to 
these  appendices  and  Figure  6 through  10  may  prove  helpful 
in  reading  the  following  subroutine  developments. 

Efferent  Subroutines . The  efferent  elements  are  com- 
prised of  all  the  subroutines  under  and  including  STRTAN 
(Start  Analysis).  ENTPARM  (Enter  Parameters)  is  simply  used 
one  time  at  the  inital  start  of  the  program  to  clear  and 
preset  the  appropriate  flags,  and  to  enter  all  of  the 
necessary  constants.  WAIT  (idle  mode)  would  normally  be  a 
constant  test  loop  for  interrupts,  but  in  this  program, 
which  was  developed  for  a multiuser  computer  system,  flags 
are  set  so  that  it  sequentially  selects  DATINT  (Data 
Interrupt)  and  CONINT  (Control  Interrupt).  DATINT  woxfLd 
normally  be  an  interrupt  routine  to  read  in  a single  point 
of  xy  information,  but  for  this  program,  it  is  used  as  a 
simulation  routine  which  generates  a full  input  buffer  (120 
points)  of  simulated  received  signal  space  coordinates. 
CONINT  would  also  be  an  interrupt  routine  to  read  in  ex- 
ternal display  selection  requests.  In  an  expanded  system, 
the  program  would  operate  continuously  and  process  all  data 
as  the  data  is  received,  and  automatically  print  out  only  a 
very  limited  amount  of  processed  data  according  to  in- 
structions which  would  be  set  up  within  the  RECRDAT  (Record 
'Data)  routines.  CONINT  would  be  used  to  request  the 
additional  current  data  or  summary  printouts.  Since  a 
current  data  analysis  system  does  not  require  continuous 
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operation  of  the  program,  CONINT  would  be  set  up  as  a key- 
board servicing  routine  which  would  allow  an  operator  to  set 
the  appropriate  display  request  flags  and  the  program  start 
flags.  But  in  this  program,  CONINT  is  used  as  a routine 
which  internally  sets  the  flags  to  have  the  current  data 
displays  printed  out.  DATRDY  (Data  Ready)  is  a subroutine 
that  transfers  an  input  buffer,  when  it  is  filled  with  120 
points,  into  another  memory  location  where  the  xy  data  can 
be  processed  by  other  subroutines  without  affecting  the 
input  buffer.  In  the  expanded  system,  this  subroutine  would 
be  quite  critical  and  timing  would  have  to  be  carefully 
analyzed  to  insure  that  data  would  not  be  lost.  However  in 
the  current  data  analysis  system,  the  subsequent  data  points 
after  the  initial  filling  of  the  input  buffer  are  not  impor- 
tant and  therefore  the  DATRDY  timing  is  not  critical. 

Transformation  Subroutines . The  transformation  rou- 
tines include  CURRDAT  (Current  Data),  RECRDAT  (Record  Data) 
and  all  of  the  routines  under  them.  XYTCFOL  (XY  to  Polar) 
transforms  120  points  of  xy  cartesian  coordinates  into  120 
points  of  polar  coordinates  (amplitude  and  phase).  APDEY 
(Amplitude  and  Phase  Deviation)  determines  the  amplitude  and 
phase  deviations  of  each  received  signal  space  point  from 
its  respective  reference  or  target  point  in  the  signal 
space.  FRQSPEC  (Frequency  Spectrum)  determines  the  fre- 
quencies of  the  phase  deviations.  It  is  a FFT  (fast  Fourier 
transform)  program  (Ref.  4)  which  processes  120  points  with 
a 128  point  FFT.  The  8 additional  points  in  the  transform 
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are  filled  with  zeroes.  COMPRES  (Data  Compression)  and 
HISTORY  (Permanent  Data  Storage)  are  routines  intended  only 
for  future  expansion  purposes.  The  COMPRES  routines  would 
consist  of  complex  statistics  routines  to  compress  the  data 
by  determining  averages  and  standard  deviations  of  amplitude 
and  phase  deviations,  phase  jitter,  phase  dispersion,  noise, 
number  of  phase  and  amplitude  hits,  and  line  quality.  The 
HISTORY  routine  would  prepare  the  compressed  data  for  perma- 
nent storage  with  time  and  line  references  for  cataloging 
purposes,  then  store  the  cataloged  data  in  permanent  files. 
These  two  routines  were  not  developed  in  this  thesis  as  they 
are  not  required  for  the  current  data  analysis  and  display. 

Afferent  Subroutines . The  afferent  elements  are  com- 
prised of  the  subroutines  under  and  including  DISPLAY. 

SELDAT  (Select  Data)  and  SELDIS  (Select  Display)  are  selec- 
tion routines  which  determine  through  the  status  of  control 
flags  the  data  routines  and  display  routines  which  are  to  be 
executed.  The  SLCURR  (Select  Current  Data)  retrieves  all  of 
the  current  data  required  to  plot  the  three  current  data 
displays  of  signal  space,  amplitude  and  phase  deviation,  and 
phase  deviation  frequency  spectrum  and  it  puts  the  data  into 
a file  which  can  be  used  by  each  of  the  current  data  display 
routines.  The  3LFMSM  (Select  Five  Minute  Summary  Data), 
SLFLN3M  (Select  Line  Summary  Data),  3LHRSM  (Select  Hour 
Summary  Data) , and  3LDASM  (Select  Day  Summary  Data)  routines 
are  only  routines  included  for  expansion  if  the  RECRDAT 
routines  are  developed.  These  four  selection  routines  would 
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obviously  have  to  fit  in  with  the  cataloging  procedures  in 
the  development  of  the  HISTORY  routine.  The  three  routines, 
DISXY  (Display  Signal  Space),  DISAPD  (Display  Amplitude  and 
Phase  Deviations),  and  DISPDFS  (Display  Phase  Deviation 
Frequency  Spectrum) , are  the  line  printer  routines  which  put 
the  appropriate  current  data  retrieved  by  SLCURR  into  a 
formatted  output,  then  direct  a printout  of  the  current  data 
displays.  The  remaining  five  summary  display  routines  are 
again  only  intended  for  use  with  system  expansion  of  RECRDAT 
routines.  They  would  format  the  retrieved  HISTORY  data  and 
direct  a line  printer  or  video  output. 

Operations  Concepts . The  FORTRAN  program  in  Appendix  C 
and  an  actual  working  prograjn  will  be  somewhat  different. 

The  FORTRAN  program  was  designed  for  use  as  a simulation 
program  on  a multiuser  computer  system.  Therefore,  it 
starts,  produces  its  own  data,  then  executes  the  analysis 
routines.  It  starts  with  TXLNAN , initializes  with  ENTPARM, 
generates  one  input  buffer  of  120  points  of  xy  data  in  the 
DATINT  subroutine,  and  transfers  the  data  in  the  DATRDY  sub- 
routine. Then  the  program  processes  the  120  points  in  the 
CURRDAT,  XYTOPOL,  APDEV , and  FRQSPEC  subroutines,  and  returns 
to  CONINT.  Here  it  reads  in  the  display  requests  (sets 
appropriate  flags),  then  performs  DISPLAY,  SELDAT,  SELDIS , 
DISXY,  DISPDFS,  and  finally  stops.  In  a real  system,  the 
program  would  read  in  a data  point  in  the  DATINT  interrupt 
routine,  then  sit  at  WAIT  until  another  interrupt  for  DATINT 
was  received.  When  120  points  would  be  received,  the  program 
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would  continue  as  in  the  previous  description.  If  the 
* system  would  ever  be  developed  into  the  expanded  system, 

then  the  program  would  have  to  run  continuously.  Instead  of 
stopping,  it  would  just  jump  back  to  the  WAIT  mode  and  begin 
the  cycle  again.  However,  while  all  the  processing  of  the 
120  points  was  being  accomplished,  the  DATINT  interrupt 
routine  would  still  be  bringing  in  new  data  points  so  all  of 
the  other  processing  would  have  to  cycle  through  and  back  to 
WAIT  within  50,000  microseconds  in  order  to  catch  the  next 
buffer  of  points.  A great  deal  of  timing  considerations  and 
advanced  programming  techniques  would  be  required  to  success- 
fully expand  the  system  into  an  actual  working  system. 
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IV.  Discussion  of  Results 

Though  the  entire  analysis  system  was  not  constructed 
and  the  analysis  program  was  only  simulated,  there  are 
sufficient  results  from  the  simulation  to  show  that  the 
current  data  analysis  system  will  provide  quantitative  in- 
formation on  transmission  line  perturbations.  The  three 
current  data  displays  of  signal  space  pattern,  amplitude  and 
phase  deviations,  and  phase  deviation  frequency  spectrum  can 
be  used  to  observe  and  measure  several  of  the  line  parame- 
ters for  perturbations  which  cannot  be  measured  on  the 
oscilloscope  eye  pattern  display.  Though  an  oscilloscope 
with  a "freeze"  or  memory  capability  can  be  used  to  obtain 
displays  similar  to  the  computer  current  data  signal  space 
display,  it  still  will  not  have  the  precision  of  the  com- 
puter display.  The  average  values  of  deviations,  noise, 
phase  dispersion,  number  of  hits,  and  phase  jitter  are  not 
directly  available  through  the  three  current  data  displays, 
but  the  exact  instantaneous  values  of  deviations  and  phase 
jitter  are  directly  shown  on  the  displays.  These  displays 
can  be  used  to  roughly  see  the  average  values  of  deviations 
and  jitter,  and  they  can  also  be  used  to  determine  noise 
levels  if  reference  displays  are  available  for  comparison. 

'simulation  of  Line  Perturbations 

The  generation  of  simulated  data  was  accomplished 
within  the  DATINT  subroutine.  The  procedure  consisted  of  a 
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DO  loop  which  cycled  120  times  to  fill  the  input  buffer  with 
120  xy  signal  space  points  as  if  they  had  been  serially 
received.  Within  the  DO  loop,  there  was  a uniform  random 
selection  of  one  of  the  l6  target  points  (which  were  con- 
tained in  memory  in  both  polar  and  cartesian  coordinates), 
and  then  amplitude  and  phase  deviations  were  added  to  or 
subtracted  from  the  coordinates  of  the  target  point  to 
obtain  a simulated  received  data  point  in  polar  coordinates. 
Then,  while  still  in  the  DO  loop,  the  coordinates  are 
switched  to  cartesian  integer  values  between  plus  or  minus 
128  for  entry  into  the  input  buffer.  The  two  equations  used 
to  generate  the  amplitude  and  phase  coordinates  in  the  DO 
loop  prior  to  conversion  to  cartesian  coordinates  were: 

A = TARGET*REAL( APT(IT))+AD+AH  (l) 

P = AIMAG(APT(IT) )+PJAMP*C0S( 2*PI*I/PJER)+PD*NPD+PH  (2) 

where : 

A = amplitude  of  simulated  received  data  points 

P = phase  of  simulated  data  point 

TARGET  = 0 or  1 multiplier  depending  on  whether  or  not 
a lineout  is  desired 

REAL( APT( IT ) ) = a uniform  random  target  point  amplitude 
value 

AIMAG( APT{ IT ) ) = the  same  uniform  random  target  point 

phase  value 

PJAMP  = phase  jitter  amplitude 

PJER  = phase  jitter  period  ( I = DO  loop  iteration) 

AD  = amplitude  deviation  (normal  gaussian  random 
variable  with  a standard  deviation  of  3IGMA) 
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PD  = phase  deviation  (normal  gaussian  random  variable 
with  a standard  deviation  of  SIGMA) 

NPD  = 0 or  1 multiplier  to  include  or  delete  random 
phase  deviations 

AH  = constant  amplitude  hit  deviation 
PH  = constant  phase  hit  deviation 
By  appropriately  changing  the  values  of  the  above  constants, 
the  line  perturbations  of  lineout,  normal  gaussian  noise, 
normal  gaussian  noise  on  the  amplitude  with  pure  jitter  on 
the  phase,  ijormal  gaussian  noise  on  the  amplitude  and  phase 
with  phase  jitter  added,  and  phase  and  amplitude  hits  can  be 
simulated . 

Quantitative  Value  of  Current  Data  Disnlavs 

Different  sets  of  simulation  variables  were  used  to 
generate  a set  of  displays  for  each  of  the  six  line  pertur- 
bations previously  described.  Comparison  of  the  .computer 
current  data  signal  space  displays  with  oscilloscope  displays 
shown  in  Figure  4 correspond  quite  well.  Also,  comparisons 
between  the  computer  displays  themselves  show  that  several 
different  quantitative  values  of  the  line  perturbations  are 
readily  available . 
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Transmission  Lineout.  The  lineout  condition  was  simu- 
lated within  the  DATINT  subroutine  by  using  data  values: 

SIGMA  =3.0 
PJPER  =1.0 
PJAMP  =0.0 
TARGET  =0.0 
NPD  =1.0 
AH  = 0.0 
PH  = 0.0 

These  values  cause  Equation  1 to  provide  Normal  Gaussian 
random  variables  with  a standard  deviation  of  3-0  for  a 
received  amplitude  and  cause  Equation  2 to  provide  uniformJLy 
distributed  phase  angles  between  plus  or  minus  180  degrees. 

The  resulting  signal  space  pattern  display  in  Figure 
llA  corresponds  directly  with  the  oscilloscope  display  shown 
in  Figure  4.  The  "T"  marks  shown  in  the  computer  display 
mark  where  the  reference  or  target  signal  space  coordinate 
points  are  located  and  the  marks  denote  the  actual 

received  coordinate  points.  As  expected,  they  are  all 
located  near  the  origin. 

Since  the  amplitude  and  phase  deviations  are  immaterial 
during  a lineout,  the  FORTRAN  program  is  set  to  identify  any 
received  amplitude  of  less  than  12  as  a lineout  condition. 
The  amplitude  and  phase  deviation  display  in  Figure  IIB 
'reflects  in  the  lineout  condition. 

The  phase  deviation  frequency  spectrum  (Fig.  IIC)  is 
just  the  absolute  value  of  a Sine  function  since  the  phase 
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deviations  appear  as  a pulse  due  to  their  uniform  random 
distribution. 
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Figure  llA,  Current  Data  Display:  Lineout  Condition 
(Signal  Space  Pattern) 
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Figure  IIB.  Current  Data  Display:  Lineout  Condition 
(Amplitude  and  Phase  Deviations) 
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Figure  IIC.  Current  Data  Display:  Lineout  Condition 
(Phase  Deviation  Frequency  Spectrum) 
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Normal  Gaussian  Noise  on  Amplitude  and  Phase . This 
noise  condition  was  simulated  within  the  DATINT  subroutine 
by  using  data  values: 

SIGMA  =3.0 
PJPER  =1.0 
PJAMP  =0.0 
TARGET  =1.0 
NPD  =1.0 
AH  = 0.0 
PH  = 0.0 

These  values  cause  Equations  1 and  2 to  provide  amplitude 
and  phase  values  which  are  uniformly  distributed  between  the 
16  target  points  with  normal  gaussian  distributions  of 
amplitude  and  phase  deviations  about  each  respective  target 
point  with  a variance  of  3.0. 

The  signal  space  display  in  Figure  12A  corresponds 
closely  with  the  oscilloscope  display  in  Figure  4 except 
that  the  received  computer  display  points  are  more  distinct 
and  the  reference  target  points  ( "T"  points)  are  provided 
for  more  clarity. 

The  amplitude  and  phase  deviations  display  in  Figure 
12B  definitely  shows  an  observer  more  about  the  nature  of 
the  noise  and  just  the  signal  space  display  as  the  average 
values  and  the  variances  of  the  deviations  are  roughly 
available  just  through  observation.  If  different  displays 
for  various  levels  of  known  noise  can  be  produced,  then 
these  could  be  used  as  standards  for  comparison  and  actual 
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values  of  noise  could  be  determined  for  operational  trans- 
mission lines. 

The  phase  deviation  frequency  in  Figure  12C  shows  the 
frequency  distribution  is  randomly  distributed  between  0 and 
1200  Hz  with  no  dominant  frequencies. 
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Figure  12A. 


Current  Data  Display:  Normal  Gaussian 
Noise  Condition  (Signal  Space  Pattern) 
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Figure  12B.  Current  Data  Display:  Normal  Gaussian 
Noise  Condition  (Amplitude  and  Phase 
Deviations) 
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Figure  12C.  Current  Data  Display:  Normal  Gaussian 

Noise  Condition  (Phase  Deviation  Frequency 
Spectrum) 
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N ormal  Gaussian  Noise  on  the  Amplitude ; Pure  Jitter  on 
the  Phase . The  noisy  amplitude  and  pure  phase  jitter  condi- 
tions were  simulated  within  the  DATINT  subroutine  by  using 
data  values: 

SIGMA  = 3.0 

PJPER  =24.0 
PJAMP  =10.0 
TARGET  = 1.0 

NPD  = 0.0 

AH  = 0.0 

PH  = 0.0 

These  values  cause  Equation  1 to  provide  amplitudes  with 
noise  as  in  the  previous  case,  but  they  cause  Equation  2 to 
provide  phase  angles  with  a pure  phase  jitter  of  amplitude 
equal  to  10  and  a period  of  24  received  points  or  100  Hz  at 
2400  baud. 

The  computer  signal  space  pattern  display  in  Figure  13A 
shows  roughly  that  the  phase  deviations  are  arced  around  the 
target  points  as  in  the  oscilloscope  display  in  Figure  4. 
However,  as  noted  earlier,  addition  of  some  noise  on  the 
amplitude  spreads  received  distribution  around  the  target 
point  in  such  a way  that  phase  jitter  is  difficult  to  see, 
much  less  evaluate. 

The  amplitude  and  phase  deviation  display  in  Figure  13B 
clearly  shows  the  phase  jitter  and  phase  deviation  frequency 
spectrum  in  Figure  13C  shows  the  principle  frequency  at  100 
Hz . 
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Figure  13A.  Current  Data  Display:  Normal  Gaussian 
Noise  on  the  Amplitude  with  Pure  Phase 
Jitter  (Signal  Space  Pattern) 
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Figure  13B.  Current  Data  Display:  Normal  Gaussian 
Noise  on  the  Amplitude  v/ith  Pure  Phase 
Jitter  (Amplitude  and  Phase  Deviations) 
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Figure  1 3C . Current  Data  Display:  Normal  Gaussian 
Noise  on  the  Amplitude  and  Pure  Phase 
Jitter  (Phase  Deviation  Frequency  Spectrum) 
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Normal  Gaussian  Noise  on  the  Amplitude  and  Phase ; Plus 
Phase  Jitter.  The  noise  and  jitter  conditions  were  simu- 
lated within  the  DATINT  subroutine  by  using  data  values: 

SIGMA  = 3.0 

PJPER  = 3.0 

PJAMP  =10.0 
TARGET  = 1.0 

NPD  = 1.0 

AH  = 0.0 

PH  = 0.0 

These  values  cause  Equations  1 and  2 to  provide  amplitude 
and  phase  values  with  noise  as  in  the  previous  example,  but 
Equation  2 also  adds  a cosine  jitter  component  with  magni- 
tude 10  and  period  of  3 received  points  or  800  Hz  at  2400 
baud. 

The  computer  signal  space  pattern  display  in  Figure  14A 
hardly  shows  the  presence  of  phase  jitter  and  even  the 
amplitude  and  phase  deviations  display  in  Figure  14B  fails 
to  obviously  show  phase  jitter  (though  it  does  show  some 
consistent  deviations  at  each  extreme).  However,  the  phase 
deviation  frequency  spectrum  in  Figure  14C  shows  a large 
frequency  magnitude  at  800  Hz  which  would  be  faster  than 
the  human  eye  could  even  detect  on  the  oscilloscope  display 
of  the  signal  space  pattern. 
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Figure  14a.  Current  Data  Display:  Normal  Gaussian 
Noise  on  the  Amplitude  and  Phase  with 
Phase  Jitter  Added  (Signal  Space  Pattern) 
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Figure  14B.  Current  Data  Display:  Normal  Gaussian 
Noise  on  the  Amplitude  and  Phase  with 
Phase  Jitter  Added  (Amplitude  and  Phase 
Deviations) 
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Figure  14C.  Current  Data  Display:  Normal  Gaussian 
Noise  on  the  Amplitude  and  Phase  with 
Phase  Jitter  Added  (Phase  Deviation 
Frequency  Spectrum) 
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Phase  and  Amplitude  Hits . The  phase  and  amplitude  hits 
were  simulated  within  the  DATINT  subroutine  using  data 
values : 


Amplitude  Hit 

Phase  Hit 

SIGMA  = 

3.0 

3.0 

PJPER  = 

1.0 

1.0 

PJAMP  = 

0.0 

0.0 

TARGET  = 

1.0 

1.0 

NPD  = 

1.0 

1.0 

AH  = 

10.0 

0.0 

PH  = 

0.0 

10.0 

These  values  caused  Equations  1 and  2 to  provide  noisy 
values  for  amplitude  and  phase,  and  they  also  added  a 
constant  deviation  to  the  ajnplitude  or  phase  to  simulate 
the  hit. 

The  computer  signal  space  pattern  display  in  Figure  15A 
and  Figure  l6A  respectively  show  amplitude  and  phase  hits  as 
in  the  oscilloscope  signal  space  patterns  in  Figure  4. 

The  amplitude  and  phase  deviation  displays  in  Figures 
15B  and  16B  show  the  quantitative  values  of  the  hit  devia- 
tions. They  also  provide  a rough  indication  of  the  mean 
deviation  (lO.O  in  both  sample  displays  of  amplitude  and 
phase  hits) . 

The  phase  deviation  frequency  spectrum  simply  shows 
'that  the  phase  deviations  are  varying  randomly. 
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Figure  15A.  Current  Data  Displays:  Normal  Gaussian 
Noise  Plus  an  Amplitude  Hit  (Signal  Space 
Pattern) 
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Figure  15B 


Current  Data  Display:  Normal  Gaussian 
Noise  Plus  an  Amplitude  Hit  (Amplitude 
and  Phase  Deviations) 
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Figure  15C.  Current  Data  Display:  Normal  Gaussian 
Noise  Plus  an  Amplitude  Kit  (Phase 
Deviation  Frequency  Spectrum) 
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Figure  16A.  Current  Data  Display:  Normal  Gaussian 
Noise  Plus  a Phase  Hit  (Signal  Space 
Pattern) 
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Figure  16B.  Current  Data  Display:  Normal  Gaussian 
Noise  Plus  a Phase  Hit  (Amplitude  and 
Phase  Deviations) 


6? 


cu^jr>iT  •«?*  •MUf  '•fourvc'T  i.T'i*'*  i » ''•th  it>t^/TT  I Tt^ri  ifi.tJ.vt. 


r»r  Tir»i-T 
•i'.-n’in 


ftv*rmr.t  i-*  m •»? 


1 

1 

1 

t 

1 

! 

1 

1 

1 « 

*- it 

• 17tl 

* 

t 

1 



i 

j 

1 

1 

1 

1 

f 

t 

• 

• i. 

• « *1 

' 

I 

• t » • t • 

1 • • • 

C m •flj  «Bt  <9«  ItfOQ  1*01 


Figure  16C.  Current  Data  Display:  Normal  Gaussian 
Noise  Plus  a Phase  Hit  (Phase  Deviation 
Frequency  Spectrum) 
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Limitation  of  Current  Data  Display 

Though  the  current  data  displays  offer  more  qualitative 
information  than  just  the  oscilloscope  display  of  the  signal 
space  pattern,  there  are  still  several  limitations  and  im- 
practicalities  in  the  use  of  the  current  data  displays.  To 
determine  noise  levels  and  relative  transmission  line  qual- 
ities, comparison  sets  of  displays  with  known  values  of  per- 
turbations are  required.  Also,  the  amount  of  paper  and  time 
utilized  in  printing  the  displays  prevents  them  from  being 
utilized  constantly  as  a monitoring  system,  and  consequently 
line  perturbations  which  only  happen  sporadically  will  be 
difficult  to  display.  Therefore,  the  current  data  displays 
are  not  sufficient  for  detection  of  perturbation  trends  or 
for  detection  of  occasional  problems.  The  current  displays 
are  best  suited  for  observing  transmission  line  perturb^j.- 
tions  when  a problem  is  already  known  to  exist  on  the  line 
and  it  needs  to  be  identified  and  measured. 
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V . Recommendations 


Though  the  current  data  analysis  system  developed  in 
this  thesis  has  some  limitations,  it  would  provide  real  time 
quantitative  information  about  transmission  line  perturba- 
tions that  could  not  be  obtained  from  other  available 
monitor  systems.  Therefore  it  would  be  worthwhile  to 
implement  this  system,  especially  since  all  of  the  necessary 
equipment  would  already  be  available  within  the  RADC  labora- 
tories. Also,  it  would  be  easier  to  initially  derive  the 
minicomputer  program  from  the  FORTRAN  program  and  develop 
the  keyboard  controlling  routines  with  just  the  current  data 
modules.  Although  the  program  structure  was  developed  with 
an  expansion  capability  included  in  the  basic  design, 
several  considerations  that  were  not  critical  in  the  current 
data  system  would  have  to  be  considered  in  programming  the 
expanded  system.  If  the  current  data  system  was  imple- 
mented first,  then  the  extra  considerations  for  the  expanded 
system  could  be  defined  more  clearly  at  the  beginning  of  the 
expansion  design  and  development. 

The  design  and  development  of  an  expanded  system  would 
first  require  determination  of  all  the  operations  which 
would  be  required  to  convert  the  data  from  the  current  data 
modules  to  appropriate  values  for  a permanent  history.  Then 
the  program  structure  to  accomplish  those  operations  would 
have  to  be  developed  along  with  considering  all  of  the- 
additional  timing  constraints  which  would  arise  from  the 
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additional  processing  operations.  Just  the  determination  of 
the  operations  required  to  provide  values  for  means  and 
variances  of  aimplitude  deviation,  phase  deviation,  phase 
dispersion,  noise,  line  quality  and  to  provide  the  number  of 
severe  or  medium  phase  or  amplitude  hits  would  be  a complex, 
time  consuming  task.  Then,  to  figure  out  how  to  get  it  all 
done  in  only  49,920  microseconds  (the  time  frame  between 
input  buffers  of  120  points  each)  would  be  another  monu- 
mental piece  of  work.  Also,  the  cataloging,  filing,  and 
retrieving  routines,  plus  the  display  routines  would  have  to 
be  worked  out  for  the  different  summaries.  The  IBM-LQM 
system  was  developed  with  the  expanded  capabilities,  but  it 
was  developed  by  a team  of  engineers,  statisticians,  and 
programmers  over  a seven  month  period  (Ref.  l).  Though  such 
a system  could  conceivably  require  several  thesis  efforts, 
the  results  would  be  well  worthwhile. 
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APPENDIX  A 

Structured  Program  Development 


Cons  tunts 


i^ure  A-2.  Third  Level  Efferent  Block 


igure  A-3.  Third  Level  Transitional  Block  (Current  Processing) 
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i^ure  A-4,  Third  Level  Transitional  Block  (Permanent  Record) 
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Figure  A-5.  Third  Level  Afferent  Block 


APPENDIX  B 

FORTRAN  Program  Flow  Diagrams 
(Calling  Routines  Only) 


”70 


1 


Figure  B-5.  Subroutine:  RCRDAT 
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XV^  (•  » =C^°LX  ( ♦2‘*.  !»?<♦.  ' 

XVT(i  ) iCMPLX  (tZI.. 

XVT (' ) sCvoLX ( -24. ,-3U.} 

XVT  ( IC^•PLX  (-2'*. 

XV’'  (■>)  =C10LX  (*4^  . , H.«  .) 
*'^T(lP)=riOLX(*4i.,-u'<,) 

'^’(11) =ri”LX( . ,n*8.) 

XVT  (<  ->)  =C*lOLX<  . ,-1.  8,  ) 

XVT  (i-f ) -r:-i.-'LX  r 49  j .,C  . ) 

XVT  (1 1.)  =r'4PL  X(  H.1  0 . ) 

XVT  (1  r)  SCIPI.X  ( -m  . ,5  . ) 

XV’  (t  fi)=C'1PLX  (-9j  .,0  .) 

C CnNV''5T  Tji^r.rT  xv  T9  TftPT.FT  A®  COn'93INATe«; 
10  1 T=l,l& 

A = r.*  0$  {xyt  (I) ) 

Vr’i-AL  (XVT  (I)  ) 

v=a*M^(;{XY’'(i)> 

° = '•9!  7’95*flT(fj5(  Y,x» 

SPT(T)=CVP1.X(A,>) 

I ■n'r’N'iE 

r fjrr.ATP  fiTPfi  oA9A*i!CT"R  plag 
r='HT">aP  = N0 


"TJT’N 

.r:4'y 


lut 
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MO  H*  O O O 


SlinooilTINE  HftTT 


WAIT  FO’  INTER'^UPTS  (SIMULATfq  /\fit) 

•JOMOLFy  XVT  (If),  ixy  < ,XY  fl?C  ,2)  , PFXYT  (16)  ,PFXY(12  0) 

:0MOl  FX  AOT(16),lP(12J,2),4P0(12D,2),P0FS(120,2) ,PFAOO(120) 
''OHPLFY  PFP'1FS(1?0) 

IMTr-FD  LN.LN'-FL,  YFc:,  ►n.LNC 

T.  tf-.F"  FrONTtJT,Fr)ATTf|T(2)  , "i  1 ( ’)  , FL I ( 2 ) , FOT  ( 2 ) 

-•■Tr'-,F7  F-TAIFT  Y,fp.SEL.FINI«:«,-=''(TPAP,FIXYFIIL(2) 

T.)Tr-rr>  FOLXY.Fp.  pIFS  , F0LF«? m, FPl LNS H , FPLHPS-^,  FPLO A<;m , F^L  A PO 
P'^tL  'I'fLNTI'iF  (?) 

'“OMMON  XYT,IXY,<Y,PFxyT,PFXY 
roiPTU  A°T,AF,A’0,pnFs,pc-apn 
PFPPFS 

LN.LN'^EL,  YPP.-ITtUH" 

:T^M'1I)  Ff^n'tlMT.FDATItlT.FlO.FUl  , - DT 
“OM-nn  FPATorv,  - P.PPL  ,FT-lI';‘',pp'JTP4P.  .FTXYFUL 
:T“‘1ph  FPLXY,F=Lp"'F';  ,FPLPmp“,FplLNSH,FPLm^SM,FPLOASM,FPLAPO 
'D-'Mifj  nTtLfiTii-:  ,iayt,tymf 
n'l-IM'IM  PL0TMAP(  I ID,  I’lO) 

r)rTFP‘*T*(P  OATA  P-  ^JN^PPL  IMTFpP')°T 
TP(FPATINT{1)  .ET.YFS)  '.O  TP  2 
IFfFTATINT(?)  .CD.  ypp)  '-.n  Tf)  ? 

TF(F''nK|If)T  .FO.YFS)  '■.n  ’O  6 

■'P  1 

prA-  Tig  "4TA  POINT  XY(LfJ) 

L-t  = l 

30  TP  i. 

T .LN='> 

<»  OALL  OATIMT 

TP(F:xyFUL(UN) .EO.YFS)  GO  TO  6 
GO  TP  1 ■ 

0 PFAT  IIJ  roNTPOL  (OIsPlAY  GFLFOTTON  op  parameter  GRANGE) 

"Alt  rONINT 

6 3-*. pit 

rip 
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cj  e 


SUnpoOTlNP  CONINT 


PROCFIS?  PONTROL  INr-TRRUPT  (5IMIJL4TI0N  AT  AFIT) 

;OMPlPX  XYT(16) ,I XY( 1?3 ,3) ,xv (i?o ,2) , prxVT (16) ,Prxv(l?0» 

“O'lOi.FX  4PT(1‘), 10(130,  ?),Aon(l?3, 3)  , POF'l  ( 1 20 , 3 ) ,PF&PP(130) 

'■OMPI.  FX  PFPDF''(13:) 
tsTr^ro  LN,LNFFL,Y-'',NO,LMf' 

I'lTF-F’’  FCOltIMT,  (3),FtO(2),FLl(?),FnT(3) 

T 'T""FO  FTflTR0Y,FD|.SFL,FIMT3M,FF>jToiP,FIXYrUL(3) 

T-'rr-eo  FPLXV,Fn_  pipe  ,rrLr‘ie',.,roLLf43M,FPLHOS‘1,FPLOA<;M,Fi'L4PO 
ofil  0T,LNTImF(3) 

■'  I'l'tOM  XYT,  IXY,  < Y,  oFXYT,  OFXY 
APT  ,40,400, P''FS,OFiiDn 
r.OM'I'N  PFPOFF 
'iTMMON  LN,LN' FL,  YFO  , HT,LUC 
n IMOH  FC0M*'T,FO1TInT,FL0,FuI,f0T 
-OMM  IM  F01T40Y,  FP1.PFL  .FT. |Tc:,4,='r,jjP4C>,  pixyful 

pplxy,  fp.  3-!F«;,rp|_irH3f',FPLLMSM  , FPL HSSH , FP  LO ASH  , F'-LAPO 
''T-I'I'IM  nr  ,L(:TIH-' , Uv  T,Tv  •<F 

TOHHou  plothap( 13:, 13:) 

r;  <;p*  npp'pFO  TISPLAY  FLAFS  AUO  O^pT^po  L I HF 
F'>_yv  = YES 
FPLflnnrYfS 
coi_o-ir^sYFS 

LUFF!  rl 

C «TFT  Aro  popIATF  CONTROL  FLAGS 

. I.'  (FPL  XY*F'»LAP0»FPLPOFS*FPLF“'3H*FPLLNSH»fplmps,1*PPL0ASH.GE.YFS)  go 
l^’O  1 

.fplS-L-NO 

TO  "0  3 

1 fplS-'l=YFS 

2 Fr-nutUTsUO 
fT»JT''H  = YFS 
OF-'I-’U 

r-jo 
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CJ  tj  tJ 


<;i)OpTiniNE  OATINf 


INT^pail'^T  POUTINF  TT  ?EflO  TI4  XV  1ATA  POH4T 
(TFMnOP.’PY  TEST  ROUTINE) 

SOMPi.rx  XYT(lfj),IXr<l2i?,?),XV{l?C,?),“FXTT(16),PFXY(l?0) 

SO'ini  FX  A'’T(l^),^P(l?a,2)  ,APn(120,2)  ,POFS(  120,2)  ,PFAPn(t20) 
:n“'’LFX  PFPOFS(l?a) 

T'lTrrr?  LN , IN'"  EL  , Y^’S  , HO 

TMrrFo  FCOHIMT,Fr)!VTTT  (2)  ,FL0(2),FLI(2),fdT(2) 

: ,TFr,rr>  FO S T Pn Y,  f 5; tl , FTHTSH  , EN TPA  R , FI < y Fiju  ( 3) 

TMTF-rn  FPLXY.FP. “OFF , ffl Fm^ m , fol lnS M , FPU E S M, FPLOASM , FPLA PO 
p"Al  '’’•^lhti-if  (3) 

XYT,  IXY,  < Y.PFXyT.ofxy 
COMMON  AOT.AF'.A’ OtPOFEjOFAPT 

OOMMPfl  ppPriFF 
EO^YP.M  LN.LHFEL,  YEF,‘IO,LNC 
'.T'V'I*!  FCOUI'IT,  ^OATT'IT  .Flo,  flt  , rpT 
"OMXPH  FOATROY,  FOLCf  L,p|,jj5u,FrNTPAR  .FIXYFUL 

FPLvy.FP.POFf  ,FpLFHFM,FP;.LNSN,FPLHRS«,FPLOASM,FPLAPn 
“INMOH  OT.LNTHE  ,iAYT,''Y'1F 

--)MHrn  PLOTHfC(  I t],  ]■■;!;) 

C r,F>)FP4TF  F^fiPLE  peCETYfo  POINTS 


O'^AU  AVE,SIGMA,PJA'*o,P  )Pfp  , T A P',  £T  , NOO  , AH,  PH 


•5“  ft 

AVF 

^’an'  f 

SIGHA 

■ ''■■ftn*  f 

0 JOFR 

H'  , 

PJANP 

TAOGET 

NFO 

9rr f 

AH 

DC*  4 » 

PH. 

' T'^  120  SAMS.F  PFCfjvfp  oaTA  POINTS 

1 .1'92 

1‘  1 = 1,120 

0 Of^jfp'tf  panOOH  UNIfopH  TARF.fT  FFlfSTION 
♦lfi*RANF(0) 

C r.F^ir.-?' TP  »anOOH  SAU5SIAN  AHOLITIIOE  DEVIATION 

Aoi  - 

00  -O  J=l,12 
0')"  = =Af|F  (AO) 

AO-  A''  + ,)U'' 

’O  0NTT»||)F 

AO=  ( A-’-r,)  ♦SIOMA-  AVF 

" OFNr'”T''  PAtJOOM  GAUSSIAN  PHASF  OFylATION 
p"  = 

0-  'T  XtI,!? 

■ I ' - c A T ( r ' ) 
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Of^:  PT*nO'<H 
30  COtrT'JUF 

SsTft5r,FT»l<F4L  (»pr  ( I T ) ) 

OsBI-'flr,  <4PT(ITn*  PJfl1O*C0$(  ’’Ol^r/OJPER)  ♦O0»NPn  + PH 

Tx=fl  *no'(Pi»p/no) 

IV  = '»''IN(PI»P/U0) 

X:TX 

.■  = IV 

:r("p';(x)  .iE.!?!.)  X = D 
TfC  P*:  «y)  .GF.  .)  V=c 
Txv(  T,  i)=p‘<Pt  x(X  ,y) 

10  ■:n»iT’M'ic- 

r CL‘’4‘  ■iA’’a  I'lTEPOiI’T  FLAG 
=-nf  TTtj’'(LN)=Nn 

R GCT  ffU*  OiJFFE=>  FJL-  flag 
FTXV-UL  <LN) =YFS 


3CT,j-  t| 

r»|P( 


MVi 


LL  l,v/l 
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a n 


onpriufiNE  oATKor 


H07'’  IN®"’’  naift  SUFFER  IXY  TO  CO’.RFNT  ''4TA  T49LE  XT 

:;OMOLrx  XYT(lC,),IXr  ( 123,P)  ,XYri?0  ,?)  , PFXYT  (16)  ,PPXY(iaO) 
:nMPi.FX  4PT(la),^o(12i),3),.>lon(l2C,2)  ,PnFS(12a,2)  ,PF&Pn(120» 
~OMOL‘’X  PFOQFSdTC) 

Tn*c-r->  LN,L^lSEL,Y-';t’^'1t^'^r; 

’ ■ FCONIMT,  ‘^Oa’'IT  (?)  .'•LT'’)  ,FLI  (2)  .FOT  (2) 

F0ATRDY,'•OLS‘^L,FI*|IFH,'^E^^T^^i^,FI<YFUL(2) 

T>iTr-.  rp  FPLVY,P=.  “IF S • ■^PL'"  « , "PU  LN3m  , FP  . H FP LO  ASM  , F^L  APQ 
5"'L  n-,LNTIM"(?) 

'OXMnf!  XYT  , IVY,  < Y,  PPXYT,  PFXV 

:oMuon  aPT,ir',s=>n,POFS,PF5oo 

^d-hph  ocotf'^ 

COMMnM  LH,LHStU,  Y--F,'|T,LH: 

*o  *MT)  ffoni'jT,  "naTTijT  ,flo,  flt  , - ot 
?0(M1N  FOATPriYjPOLSFLtFI'irSH,- E'JTPAR.FIXYFUL 
nM-iv  pplXY.FP.  P1S'S,''PL- M3-,FPLL‘';'1,FPLMRS‘1,FnLOASM,FPLAPO 
lOMMOM  nT.LNTHF.TflYT.Ty-- 

PL')TMaP(  I Tu,  I'Y?) 

C full  imp'j’’  i/fffp  into  ftl''  xy 

IT  1 ■'  = 1,1’0 

XY(v,LN)=tXY(I,.N) 

1 ''0‘JTfmiif 

C r.L'^A'’  full  giJFFFR  flat, 

' XYFilL  (LN)  =NTJ 

C FF’  TA’PPY  FLAG 

, FMaT’FlViYES 

Rft'JPN 

Fk|ri 


c c c 


SUIRmiTINE  CURRO^T 


''U'RFNT  OaTA 

CO«=LFK  XYT  tl<>),r  xy(  120,’)  ,XY(1  ’0 ,2)  , PFXYT  (16)  ,PFXV(130) 

‘^X  APTClb) , ) P(120,2)  ,a=n  (133  ,2)  , ( 1 20, 3)  .PFAPO  (120) 

FOMPLFX  PFPnFSd’O) 

T •jf-cn  lN.LMSFL,  Y-';,f|'1,LNC 

’••jTC'.FR  FCCNINT,- O'TIMT  (3)  ,'=’L3'2)  ,FLI  (2),PDT(2) 
f.jcr- FP  criATRPY,  - PL^FL,  FINT'H.cf.jtP.IR  ,FTX  YFUL  ( 3) 

IU’’F' FP  FPlXY.FPl  piF';,FP(.fm'^u  jFP(_L,|<;h  ,FP..  ^P,Sf1,FPLnASM,FPLaPO 

PFflL  T'.L'flM-' (3) 

3T^’nH  XYT,IXY,XY,PF)(YT,oryy 
APT, BP, APT, POPS, PFAPO 
OOMH-'M  DrOTFF 
-TjM.F,,  L‘i,  LN'^fl  , Y-'P,  ‘lO.U'JO 
OT-'MOU  FCO'ir  T,  FOATi.ii  ,FLO,  Fir,  FOT 
OT'MON  FTATPHY,  pplSFL,  FTNrPH.FF'fTPAR.FIXYFUL 

''T‘‘)OH  FOLYY,  FP.  P)PF,FPLFMPP,FPLLM3H,  FPLHPSH,  FPLOAS"),  FPLA  PO 

•>0M-4-i,|  nT.LNTHF.TAY'T.FVMF 

'otM-M  flotm6P(  1 p; , 1 3: ) 

C FLFAr  PBTA  RFAOY  FlAO!  3fT  LMO  (lIXF  * PEIN5  PROCFS^FO) 

FOA’--ov  = r)0 
LN~  = ' M 

C I'DMYFP.’’  nAPTPSIAN  TO  POLAP  OODPOTTA  T^S 
'•ML  XYTOPOL 

C nAL'Ol'LATF  AmolITUDf  ASP  oMAPF  O'^YIATtON 
.FALL  a=nFY 

C falfuLATF  phase  deviation  f=FH|f,|Cy  SPFCTPUM 

FALL  For)FDC(; 


• |0 


r-  ^ V ' • 

h •,  ^ i 

■_J*  I _ i'i  li  A Aii  'k^ 
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o c o 


«5II'^R'I0▼I^E  XVTOOTL 


Cn^n/FOT  '•ft^TFSIAN  TO  POLAR  COOPTINSTFS 

~0^r>Lry  XVT  (1:0, 1 yT(  1’3,3>  ,yv  (1  20  ,2)  , PFXYT  (16)  ,PFXV(120) 
OOMOi.ex  APTdf  ),1P(  12  0,2)  ,APn  (120,2)  ,0  025(120,2)  ,PFAP0(12a) 
rOHPLfX  PFOOF':  (1  T 0) 

T'ITr'-.ci»  Lt.,LMFFL,  Y"':,‘JO,LNC 

V^Troro  FCONInT,  *r)l’’TN'^(3)  , '■LO(?)  ,FlT  (2)  , 00T(2) 

7 .jTrr.co  FOATRPY,FOL3FL,FIMr5M,FFNTOAP  ,oiXYFUL(?) 

I.jTr- ro  polXY,FO.  ’OFS  , OPL or?** ,pol LN5  1 , FFl  RPSO,  pol 0 A5'1  , FpL  A PO 
POAL  a'’,LNTI«F  (?) 

YYT,IYY,<Y,PC-YYT,=o<y 

COMTM  APT,AP,A?3,PnF5,PFAO" 

PFPOF" 

"0  (“Of)  LN,LN5£L,  YO5,N0,I.NC 
■;n'<'n*i  fconimT,otati;(T,flo,fl!  ,-DT 

•?))  FOATPr-y,  roLFTL  ,FT‘(1 5-','^  o^rpu!^  , FI  KYC-UL 

folyY,P“.  P iP'^ , =’PlFm~«,FPLLN5M,PPLHOSM,  PPLOA^M,  FdlAPO 

■'OIMIM  nT,LNTI  t-  ,OA  YT,TyMF 

olOT“11p(  1 ?C,  I'TO) 

'O'lY'P'TON  LOOP 
00  1 Tsl,123 
\:ri 05  (XY (I,LHC)  ) 

IFTcpi’M"  LINPO'IT/L  IMFIN  CONOTTION 
'Pd.po.o.)  r,o  ’•n  1; 

-PC' .LF. 12.0)  50  TO  11 
PLT(L'JC)=YPS 
50  *n  12 

0 n«TTNP  A = 0 VALUES  A NO  PET  FLO(LNO) 

10  1*0.  F=a. 

PLOd  *'0 ) =YPS 
50  '0  1? 

11  '■LOfLMOrYPP 
50  -o  12 

''  OHISP 

12  X=P"AL (XY ( r,LNC) ) 

v=a’’‘ug(xy  (i,lm:)) 

P*p'  . ■’96T795*AT1N2(  Y,X) 

0 OPCATO  AO  hole 

IT  AP(',LN(,)=CPPLX(A,P) 

1 ' 'ONT'HUe 

- --M- H 
r-)0 
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r>  <r% 


SU'^RTIItine  ^pr\r.\l 


r.  nFTPTHiH'’  amplituoe  Afn  phase  nr\f'^rTr)n 

■^dMPl  FX  XYT  ( IF  ) ,T  XV  ( 1 ?C  ,?»  ,xv  (1 ’c  ,?)  , »FxrT  (l«i),PFXV(l?0) 

COMPi.FX  APT<16)  .APd^O,?)  tAOPti  ->(;  ,7)  , POPS  (120,?)  ,PFACO  (I’O) 

2n‘<PLPx  pppoF^d’o) 

T-iTp'-.pf'  lniLn-el, yep, *io,L>r 

jurr-.ro  F.',CNI'(T,=-01T:-|’’(?1  ,'=-n('*l  ,‘^LI  (2)  ,-TT(2» 

• iTr-.PP  pnATP0V,FPLS'"L,FIMIPH,rc^XPftu,PIxVPIJL(2) 

'■  r P'.  r->  FPLXY  , PP.  0'f5,'’eL"MS‘','’°LLNSH,  FFLHPSHf'^PLnASH,  F°LAPn 

’•■aL  nT,LNTIHE(2l 

XYT,Ivv,<v,prxyT,Dr<v 
AST, AC  jA’A, PS’AO"' 

~iT  IHIH  pporiFP 

:0|Mnn  LN.LN'PL,  YFS.xn.L'l' 

On  v-'^M  FCOMl-lT,  - OaTlN’’,cu7,CLf  ,?ir 

'OMMTfJ  FTiTP"Y,cp..''-L,ciMl';M,-c«jT'’AR,PIXYCi)L 

PPLXY,  pc.  PIC  p , COL  c r3i_LS5H  , PPL  HRS**,  PPL  HASH,  FPL  A PO 

211HON  nr.LNTIlc , itYT.ryyc 

OniHPH  PLO''HaP(  1 30,  1 1: ) 

''-Ai.  P(1F)  ,nHT’4 

r,  'v;l"  *HRnur,H  ehti’f  xy  TiTLC 
■'n  1 ■'  = l,i’0 

r n-AcPK  COP  LO  ComiTTOH 

fcc-«-AL(Ao<I,Ln:))  .LE.12.)  PO  tD  3 

n '5"F'’-'If)E  OISTANC:  cPOH  PWCPV  TAPPET 

30  ’ J=l,l^ 

7<  I)=CAPS(XY(I,LNC)-XYT(J)) 

? POH'FKijr 

n nFTPPHT.ip  HiNi»*uH  OTSTAiirr 

TT'::A'iIr4i(7(i),n(r>),1C'(,T(y),D(5»,O(6),O(7),O(3>,D(q>,OC10»,O(l 

' PTCPM'*)"  tcAir<;T  rAP'-c’ 

J-:,ie 

rr  (n  nrj. )(£.'■'  (Jl  ) CO  TO  , 

c,  rALP;)'  ATP  AHPlITUlp  AnO  Pk'AS'"  IcytATION 
ac"'  < c,LMC)  = A P(  r ,lM''  ) - APT  ( J) 

c pcroi.'/p  ♦ O’-  PI  iHiToiirTY 

TP(‘  THAr,(AP0(  I,.»)P)  ) .LC.-IAP.)  APO(I  ,LN0I  =CM»LX  ( pc  AL  ( APO  ( I , L NC  ) t , 
j d'''!'-  (Ap-'(:,Lncn  TTF 0 . ) 

, 'T'J'TMIJP 

'.n  'o  1 

~ AT-;T'-,'|  ADp  /alhcc  rjp  lO  On»nT’'TO'( 

T )'■  1' •'.L'lOl  iC"OL<  ( n . , 3i.) 

t ^ ' ‘!'|C 

.ftp'y-.vv 

cHQ 


O Vi 


<;i)^orniTifiF  FROSOiC 


C 


nFTfSMIN''  '■RFOOCNry  of  PMlSi  OEVIATIOM  <FFT) 

tOM"!  ry  XVT(lb),  ' xr  ( ,yv(l?2  ,2)  , PFXYTtlf.)  .ot-xYdaoi 
«OMPL  FX  4Pr<m,A'>(12Q,?),AO"(l?!:,2)  ,OOFS(  123,2)  , PFftPn(l20l 
''()MP;,rY  pFO'ir'?  ( 1 ? g » 

TN''"*'’'’  LM,LN''fl,  x-c' ,*P,L‘J: 

T.(Tr-rr>  FrONItl T,  " TM  TNT  ( 2 ) , ='10  ( 2 ) , FL I ( 2)  , - DT  (2  ) 

’ Tf'fo  FTATRiY.s-OL^FL.FlNI^ J,=TyTD4p,FixYri)L{2) 

'•PF'FS  FFLXY.FO.piF'-  ,FPLFM':'',F3LLNS‘1,FPLH»S«,FPL0ASf1,F^LAP0 
-■"■SL  O’.lNTIX'^  (2) 

•31X>*1N  XYT,  Ixy,X  v,PFXYT,opXY 
'TMHTH  APT,AP,A>T,P0F<;,pFS3P 
20“‘'pfl  ppPIF*" 

-riMHIN  LN.LNStLf  Y-'C.NO.UN' 

'ON’^PH  Fr^ONINT,  foitINT.FLO,  F1.T  ,P0T 

FOSTtjoY,  r PlSfl  ,FlMIF-i  ,rr*|TPAR  ,FIXYF(IL 
2 0HH-.N  FPLXY.FP.  o''Fr,FPLFM':“,FPi.UNSN,FFLNRSM,FPLl3ASN,FPLAP0 
common  CT.LNTIH", 3“YT,TyHr 

"■PMOH  PL0TMAP(l3j,  130) 

20MP'.  P X PO  (12?  ) f ),Wf  T 
l(.li:P2oF 

C CY'  LF  TATfl  INTO  PI  A'^RAY  TO  PF  TRANSFORHFr) 
in  1 1=1, 1?0 

PO  ( T ) =c  ►'PL  X ( A I'd  r.  ( A PO  ( t , LN' ) > , 0 . ) 

1 ■ ^nxxnDF 

0 ■•rpo  o-MilNOFR  OF  30  array 

.00  2 T = l.?l,12fl 

oOCxirCMFLX  (C .tO .) 

2 00N’’’’MI(E 

^ .icDcp-j"  .ippoppoYAYr  OFVFPRAI,  0^  PO  ARRAY 

l = * 

n ^ T= 1,127 
TF<* . OF.  J)  GP  TO  la 
* = "■'()) 

FO'  I) =Pn(I) 

jp(')  =T 
13  < = F.' 

11  TF  j,  CO  TO  1? 

J=J-K 

<:</? 

•p  xp  1* 

!■»  ):J»X 

3 "ncr-NiiF 


P'-cpp  . •R.I'JSFOfaAT  ION 
' L = l,7 
F r • *1 


* ^ 

L.’im 
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LFl^LP/a 

U=CMOLX(1.,0.) 

Wsf;MPLX<COS(f>I/<^L')ftT  (LTD)  ,-';TN<PI/rLOAT(LEl))  » 

n''  - J=1,LE1 
no  f,  I = J,12d,L- 
If»=I+LEl 
Ta^'r(rp)*u 
nr)?  TO)spo(I)  -T 
P1(I) rPO<I) ♦! 

H Cn-fTNI'E 

.»rM'  M 

5 CON'TMUF 

4 tnUT’Ullf. 

C PUT  T^' •l'^rr«PM(:0  OAT#  BACK  TNTT  3U??ENT  OATA  FILE 
50  ’ 1=1, 1?0 
oorr(:,LNC)=PO<t> 

» ;ONTTM'IP 

~*I5 


f. r 
I r 

S \ }i  i U t-J  ^ 
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o c* 


SinPniiTiN?:  PECRO^T 


c 


Rprop'i  isTo  PFR<a  i'’'JT  ftl'’  tn  tim'’  co“<PR'’';sn  fop‘1 


:o**PLFy  XVT(lf>),iyy(l'*0f?),yv(l?^,2),pcvvT(16),PFXY(l?Q( 
:riHPLrx  a?T(iF),ap(i?o,?),ao''(i->j,?), pops (i?Ct2»fPF«PO(i20» 
?o>iP'.px  PFPnF?a?o) 

I.4Tr'-,rD  lN.LNFEl,  xrc  , ND , LMC 

••.TP'-.FP  FCONINT,  - nATTNTJ?)  .rt  K 2 ) , FL I ( 2)  . - IT  ( 2 ) 

: <TF-ro  FDaTRPX,PPL«:='LfFIHI'«f*PNTPaP,FIXVFUL(’) 

X>r<-'-FP  FPLXY.FP,  P1Fr,FPLP'<';P,-PLLHSH,FFLHRSM,FPuOaSM,Fui.APO 
Pt-BL  P’.lNTInr  (2) 

'OMMo»|  XYT  , IX  Y,  X Y,  OPXYT  ,OFX  Y 
(.PT  , tP,a’T,P'‘.'’S,PPaP'' 

*r»i‘nn  prPHFr 

“O'i'f'M  LfifLN'^ELf  Y=’S,NO,L‘r, 

pCO'JI-iT,  ’TIT  it  .PL^t  '"L’’  t -'^1 
FnaT9nv,FPL‘;FL.FI>(T':M,PFSTPaP,FIXYFUL 
~n'“*PlN  FPLXY.FP.  01P5  ,FPLP-<'H,POLLN>1,FP'.HPS*1,FPL0aSH,  FPtaPO 
nT,LNTT•^r,^BYT,■'Y'*P 
-T'-'N  PLaT“AP(  I TT,  riO) 

r I'o-ippF'r  data 
PALL  Fn*<PRF<? 


C irT?M'Mr  LIHE  peaty  FOP  pcpMANFMT  FILf 
TF<Fpt  (1)  .CO.YFS)  '.0  TO  1 
T -(P"''’  (2)  . PO.  YE3)  -,0  TO  1 

' t 

C POT  "o-*p3E<5SE0  OATA  into  PpRNANPNT  FILE 

1 full  mtctory 

PPTO’N 

FNO 


Vi  / I ' / V . T * 


cur 


a 
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c:  c c 


C0MP9i«! 


f^ri>-'P5*'5?  P«TA 

'O'O’I.PV  XYTClft),  r XV(  1?:,?)  ,»Y(13a  ,?)  ,PCVVT  (lA)  .orxYd’O) 

COMOl  PX  APT (It ) o(l?3,2) ,f an (1  JO ,2) ,oOF3< IZCt ’>  f^FApO  (120> 
20MOI  Fx  PFpnFSd’O) 

IMTF^PP  LN,LH'‘EL.  YP3,H1,L‘J'' 

TUTP'.PP  FCT'II'4T,p'1iTTtjT(?)  ,plT(2)  ,PLI  (3)  t'TlT(2) 

T TFr.FT  F0ATP0Y,pp!.$pL|PIMT3m,ppMToAR,pIXYFI)l(7) 

PP  FPLXY,FP.  P1P3  .FPLP^'IF-  ,fo^L.(3,^  .cp.^i^ps^t  FPLOAS'I.FPLAPO 
P-'AL  , LNTIf' (2) 

xYTfIYY,  <Y,CPXYT,PPXY 
nO'iMin  iPT,flc,i->n,rinF:,PFaon 
-•j.jtFAFj  r-fP3F3 

:T4‘<''ri  LM,LhF3L  I YP*^  fNIfI.M" 

FCn>llNT,  P DdlNT  ,FLn,  FLI  t - OT 
-n—(nn  F^aTP"Y,FPL'FL|F!*lT'^,PF'^’■o.^R,FIXY='UL 

'0  .Mnn  FPLXY  t FP.  OTrr,FPLF'^''’,'‘aLLN3M,FPl.HPSH,  FPLO  A3 1 f F“L  A PT 

''|m'^->fi  OTfLNT  f'^A  YTfTYMP 

-nMHOK  PlOTm1P(  i30f  !,Tj) 

JFTinFI 

•*n 


'■'PP-'IITINE  HI3T02Y 


r ui^rnPY 


■‘O^'Pl.  FX 
' T*<0'.  FX 

T.j'F'n 

TN’’"'.PO 


L P' 


;n>..iiFj 

" n ♦•*  nF| 
"•  n • ' • * 
' j^MPN 

- V*'<  ■'►I 

- VI 


YYTdf  ), 
APT  (!•')  F 
PROOFS (1 
L‘l,LfJrPL 
FF3*nf:T, 
FDATP^V, 
FPLXY , FP 
fL^UT-pc 
YYT, IXY, 
A-'T,  1P,1 
pronp  , 
LH,  L*'  - p. 

FFiVJl-'  r, 

F'lAT-'P'Y, 
FPLXY, Fo 
PT.lKTIA 
PLOTi;  ■>{ 


I yy ( 123 ,p) ,XY (1 ?0 F pFXYT (16) fPFXY  d »J) 
A“d?3F2),APP(l20,2),  OOF3  ( 120f7)f“p*P3(120) 
■’0) 

fy:s,nt,lnfc>) 

Pri'  TTNT(2)  F^l  0(2)  ,FlT  (2)  , POT  (2) 

FPLFPLFFI*ITSu,ri-,jTPlP,FlyYF|tL('5) 

. P0FC,FPi_r«^^^nLL‘J51,FPLMPS'1,  ^^1  3AFk,FPlA“0 
) 

XY.FFxvT.txy 

P ),P'FP,pr43n 


F Y-3 , ’nFL'P 
■ lATI'l'  ,PLO,  --I  T F - r^T 
- P.'  'L  fF^'ITS  J,--'|TPA-'  ,FT  <yFUL 

LPOFS,rp|_F«t;F,,FPLLNSM,  FPL  HPE.  FPL  OASI , FpL  A PO 
.PAYT.TYir 
( o . n , 1 


1 TOf  ’ 


m i.c« 


I 


.*>  r>  n r> 


SUBROUTINE  OISPL^Y 


ftFFF  = PNT  r-ONTROL  O0JTT“(F 

;;r)Mpi_ry  XYT  ( if, ) , I y Y ( 1 2 J , , v v ( 1 , 2 ) , PFXYT  (16)  ,PFYY(120) 

^O'-PL'-y  6PT<lt  ),)  0(130, BO  (120,2)  ,POFS<  120,’)  ,PF#po  (120) 
noMP'. py  PFonFS(i2  0) 

’ TF-FO  lu,ln‘-'fl,y!-«:,no,lnc 

’nrr'-.r-?  rc0MlMT,o  01TT‘l*  (2 ) , =•!  0 ( 2 ) , Fl.  I ( 2)  , "OT  (2 ) 
j,,Tr-r3  F04TFOY,oOLSFI.,Fl;|ISM,oE'(TOfti?  ,riy  YFUL  (2) 

T'l'oOFP  rPLyY,FP.  o)C«:,FrLF>lS“,ooLLN';M,FPLHFS‘1,PPL0ASM,F0L4P0 

RP4L  n'',LM*IMr  (2) 

rOYYON  yYT, lY Y, < Y, Pryyj, opy Y 

PO‘*«OH  lPT,AF,4'>0,opFS,oFaoo 

^O'tMOU  OFpriF<: 

LN,LN''EL,  y-S  , 

'OMM-iN  FCOMI'IT,  P OlTIfr  ,PL0,  FLI  ,001 

F0aTP''Y,=-D[.soL,FINI^M,ppyTP4O,oiXYFtlL 
-nMMPM  FOLXV,Fo.oiFS,ooLcyi^M,CDLLK|5H,FPLHFS*1,FPL0ftSM,  FPLiPO 
PTjLNTT'I' , 14y'T,TY'lo 
OOM^ON  PLO’M4P(  1 30,  nO) 

'“ILL  '•CLniT 
0"LL  OOLPIS 

307(1-  ti 
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SDliROUTINP  SEL0»r 


r 

r 3 

c ^rLF^TTOM  nr  selfct  laia  RonTiNc 

-nMDLTx  XYT(l^),T  XY( I’S,’) , KV (I’O ,2) ,PPXrT (16) ,PfXV(1301 
?OM3i.F<  43T(16),lO(l?0,2),4Pn(12C!,’)  .POPS  (120,2)  ,OF4PO(l20) 
rOMr*!  ry  pcpr)F<;  (I’fl) 
t«jTr-ro  lN.LNSFL,  ''-S.OO.LNC 

T .Tr'rr.  PC  ONIN  T,  - OS  T TMT  ( 2 ) , PL  0 ( 2 ) , FL I ( 2)  , ^ OT  (2 ) 

'•iT<-rr->  F04TP0Y,PfLPFL,FI’'JT"‘','"E'JTOSo  ,riyYF|)L(2) 

t.jtt-,  ro  POLXYjPO,  o->F  c , PPLP  mi;u  , pp^  LUSH  , FP^  ROS  “1,  P^L  OSSH  , FPL  4 PO 
P'^SL  ’’T,  LMTI“r  (2) 

CO»(Mni  XYT,!yY,<Y,PFYYT,3e-YY 
COH'IOfl  4PT,4r  ,4»0,P0FP,PPSP" 

CO"«OM  pponp-- 

-T.MTN  LM,LN'=FL,  yPP.'JO.L'l* 

'OMM orl  FCOUP  T,  -OSTI  .Pin,  PLI  .'■PT 

F04TPrY,'PLPPL,PT>II'H,PP'IT?4»,PI<YPIIL 
CO-tMON  FPL  YY  , FPI,  POF'' , "PLP^'P  “.FpllMP  H,  FPLHPSM,  FPL0AS>S.  FPL4PO 
PT.LNTIM-  , 04  YT,  TYMF 
FL0T'14P(  I 3';,  I'Y' ) 

TP(P”LYY»pplAPP*p“lPPPS.Op.yc'-)  -.o  TO  500 
Tr(r'i_’^Mpv,.pn.Y- 3 ) GO  to 
IFIP'LLNS  <.P0.Y-  i)  -.0  *0  502 
-rtpiLHP.'-i.PO.YP  j)  Gt  to  5C'' 

TP  (I- -1.  045-<  .EO  . YP3  ) GO  TO  o;:, 

■jG-'ir  n 


500 

' LL 

'■LCDOP 

’PI-' 

H 

50  1 

"Sl  L 

PLPHSM 

, PPTI)-' 

M 

53  2 

'ILL 

5LLKP5 

p-Tir 

*1 

50  Y 

51LI. 

5LHP5H 

50  L 

""LL 

" •»! 

M 

■"  * ;n 
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SUnROIlTlNE  SLCU^< 

C 

r «;PLrcT  ciopent  oata  po?  oi^plav 


po'1<»lfx  xrT(tf)i,r)cv(i?o,E»  t»'v(tP0*?»  ,ppvrT<is> ,ppxy(i?o) 
'“oyotry  APT(i»^),AP(i?a,?)  ,t°'^(i’0t?).P'iES(i2o,?)  .pfapikizo) 
•^OMOLfy  PF=nFS»i»0) 

I'JTr'-.ro  LH,LNS£L,YFStHO,L‘ir 

' '•r-co  FCONit  T,-0ATIHT(?)  , s’ Llf  2 ) t FL  I ( 2)  ,FnT(?) 

T.i»F'r->  FOATRrV,  ‘^PLSCL.FIMI'H.'-Fs'fpAP  ,fIi(YFUL(2> 

pjTf'.po  F3LXy,Fo.  PaFF,FPl.^M!;*•,l^OLLf^5'*,FPL^RSM,FPL0ASi,FPLAPD 

'’"AL  0’‘.L’JTIM' (?) 

-IHHOU  yYT,IXY,XY,PFxyTtPFyv 

A»r,AP,A’n,pnF';,o'’Aon 

PFPaFF 


LN,LN''.EL,  Y-'itNO.L'JC 
FCOMINT, - OATINT .fld, FLT.FOT 
FHATROY,  =’PLSrL,FINISH,c£STPAa,FiXYFUL 
-TmH  FPlXY,  Fo.PTC«;,rpLFM^«,FPLLM3«,FPLiRSM,FPL0ASi,FPLAP') 
“OMMOM  OT.LNTI'*-  ,1AYT,TVHC- 
PLOTMftPdSj,  n^) 


R r)CTr->MT-!F  QATF  Ann  TI^F 
PALL  PATFCnAYT) 

'AIL  TTHECTYHF) 

c ■:F(,rr»-  Apppop'^ia TF  cjppe'iT  data  an3  transfe!?  to  plot  files 

IF(F-'LyY.Fa.HO)  *.0  TO  12 
1-  T = l,16 
PFxyyf i>=XYT(i> 

10  ?OH'’’NI)F 

•PO  11  Tal,120 

pFifr  rT)=yY(i,LN;EL) 

11  CONT'HUE 

1?  TF(F'LSPn.EO.*lO)  GO  TO  14 
P'1  1’  T = l,120 
'>F'.-'n(  I)  =AOO(  IjUN^FL) 

IT  'nfjT’’H!)F 

14  IFCFPl  POFS .eO.NO)  30  TO  16 
PO  IP  T=1,1?0 
PrPPr-:(i)=soFS<T,LNFEL1 
IK  'DNTTM'iE 

16  pftiJ'N 

FHP 


f * ' ! It  n I 

■ - f " " t • • t 

if  fLiUii-- 
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«r>  r> 


0UM1V  5f:LFCT  0*TA  SUTROUTINFS  NOT  YFT  OFVELOPEO 


C 


•JlIPPOUTINE  SLFfl‘5'1 


complfy  x»T(i«i)  ,rxT(  IT'’,?)  ,)fr(i?0tP)  fPfxrr<i6i  ,pfxv(1?o> 

;0MPi.ry  APT(lr.),AP(123,?),*PO(123,?>,  POPS  (120, 2)  ,PFftPO<120) 
:OmOl''X  PFPOFO  (1?0) 

T irr-rp  LN,LN';FL.  YE'?  .NO.LNC 

TfjTr-Fp  FCONINT,- 01TI»r<?),Fn(2)  ,FLI  (2)  f '■DT(2) 

T-jtc-  ro  FflATPHY,  - 0L';'’1.f  T’TNT‘;h,'^c-NTP4?,pIX  YFIJL  (2) 

TMTFOFR  FOLXY , FP. P0P5 , FPLFmS H , nL LNS N , FPlmPSx, FPUOASH , FPLAPO 
P'^AI.  0'',LNTIF'^  (?) 

rO'IMTM  XYT,lXY,XY,PFXvr,Pryv 

"ONMON  apt ,AP,APO, PDFS, PFAon 

SOMMOM  PFPOFS 

COMMON  LH.LNSEL,  Y-F,N0,L'IC 

COMM  OH  Fr:0NINT,P0ATINT,rL3,riT,- OT 

common  FOAT)<r.Y,PPLFFL,FIHISM,'*FNTPAR,FIXYFUL 
"OMMON  FPLXY.FP. POrs,FPLPMSM,PPLLNSM,FPLMf<SM,FPLOASH,FPLAPn 
common  OT.LNTIME  ,0AyT,''VMr 
CntMON  PtOTMAP(  1 30, 1 TO 

rtjo 


StJPPOMTINF  SLLNSM 
r 


OOMPI.FX  XYT(l?.),rXY<12e,2),<T<120,2)  , PFXYT  ( If.)  , PFXY(  1 20) 
rOMPI  FV  APT(1»S),\P)120,?)  ,APO(t?C,2)  , POPS  (120, 2)  ,PFAPO(120) 
COMPi.ry  OFPOFSd’O) 

I'lTFOn  lN.LNSEL,  YFS,NO,LNC 

T,)TP-,r  ->  FCONINT,- OATIN’’(?)  , '‘I  O' 2 ) , Fu I ( 2 ) , FOT  (2 ) 
fl-r-.r':.  FOATP"'V,  'P.F-L,FrNI'-M,PFNTPAP,FIxYP)W.(2» 

- Tprpo  FOLXV.FPLPTFSiFPLFMSM.rPLlNStltFPLMPSMjFPLOASH.FPLAPO 
■CL  "T.^NTIMC (?) 

'OMMOM  xYT,lXY,YY,PFxv',oryy 

common  AP' , SP ,A»0,POPr ,PFf 30 

-OMMON  pfpqFS 

-O'.MON  LN.LKSFL,  Y-"S,M0,l.'jr 

-OM-roN  FOONIMT,  PO-ITI'i*  ,PLn,FLT  ,POT 

-OMM-,,)  FlATpoy,  '■FlOFL.FINITM.'-ENTPAR.FIXYFDL 

•OMmon  KOLXy  .FP.POFS,  rr(,FM^M,roLLNS  M,  FPL  MPSM,  FPL  hash,  F'LAPO 

-O'-iMON  t)T  ,L  NTI  •(■  ,OA  yTjTXMF 

?nMMON  PlOT«^o{ 13i,l TJ) 

PP’IION 

'NO 


r. 


SlinROIITlNE  StHRS-l 


SOHOL^X  XYTdf  ),txr(120t!»)  fX»(t?e,?>  .PPXXTdS)  ,PFXY(1?0) 
C;nN|PLFX  flPTd?.),lPd?3.2)  ,tP''d?0  t2)  , PnFS(l?0,2l  ,PFAPO(120» 

'omplfx  pfpof';(i?o) 
i»itf'-.fy  ln.lhtel, ye';,m'i,lmf 

T./TT-ro  fcoN'I'.'T.fdaTTV”  r?),PLO(?),FLT('>),Ff'T{?) 
tUTFr.rP  FTATPFY,  fplS'L,  FX^JTSp.-fxTPAP.fIxyFULCZ) 

; •YFCro  F°LXY  ,FP.  P3F<i,FOLF>lSM,F3LLMSM  ,FP.^?S‘1,FPL0ASM,  FPLAPO 
•■«L  OT.LMTTY-  (?l 

XYT.IXY.XY.PFXYT.PFXY 
:0-<M0N  APT,AP,A5n,Pn«'S,PFAPn 
PFPQFS 

tOMMOf)  LN,L^FfL,Y'‘F,MT,L'jr 
rOHMON  FONMT,  Fr)iTI>|T,FLfl,  F1.I  .^OT 

FOATFnY.F  plFfl,  FINISH, FEMTPAR,  FIX  YFUL 

FPLXY,  FP.  PTF?  jFPLFM'-fFpLLnSMjFP'.HRSMfFPLOASH.FPLAPO 
CO'l'nN  ot,lktim-*,.yayt,’vif 
-OMMTN  PLOT>ieP(  133,  130 

pfT'J^N 

riijn 


FIIORlUTIKlF  SLOAS-I 


'OHPt.rx  XYT  d<>),IXY<  12:  ,2)  ,XY  (123  ,2)  , PFXYT  d»i)  ,pfxY<1?0) 
^OIPLFx  APTdF,),APd23,2)  ,AP5(12C  ,2)  , POPS  ( 120, 2)  ,PFAP'Y(120> 
CO^oi  FX  FFP0F5  <1'’0) 

IN’T'-.rR  LN,LNiFL,  YFp,no,lK4C 

T*ITf^Fp  FCO'4IriT,FT8TrN’’<'’)  ,fl  1<2)  ,flI(2»,f')T(2) 

» iTP-.r^  FOATPrY,- P.':FL,FIfrXH,Pr,4TPAP,FIXYFUL(2) 

• TP'.i-p  F“LXY,FPup''p'',pplfmF'<,'’“.LN5M,FPLHPSY,FPL0AFH,FplAP0 

•>r«L  f"' ,L>JTImF  (31 

yYT  , IXY,  Y Y.PFXYT,  PPVY 
iPT,AP,AP0,oPF'‘,PF1P'' 

"O'M-'M  PFP^F~ 

'T.MOti  tN,L‘'-5'L  , Y"'-,'n,L‘l* 

*T"nM  FCONIM,  fT4T:')t,c-L0,FLT,FP: 

FTATk'  Y,FPLSrL,ritlTOP,-'’MTPAP.,FlXYFIIL 
3 0 fPLXY.FP.  "•'P3,FPLF“*'5'<,''PULN3M,FPLHO<5H,FPLnAx;H,FPL»Pn 

0T,LHTI1-  ,1AYT,TirMF 
- FLUT-toC  13'J,  3 33) 

-•■'ll  *1 
r»(n 
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lY*  T 


^ if 


r>  o r> 


SU'IRnUTINP  SELOI^ 


SFLrCTIOM  OF  QISPLAr  OUTPUT  ROUTI'^'' 

roMOt.rx  xYT(ic),  iXYdeOf?)  ,xv  (i?o ,?» ,prxrT(i<i)  ,pfxy(i20» 
COHPI.  rx  APT  (If  ),  JO(l?J,2)  ,»00  (t?0  ,2)  , POPS  1 120t2)  , pFAPO  (120» 
tOMPl.FX  PPP0F<:  (1?0) 

TUTf'.F^  LN.LNSEL.YFS.^jn.LN' 

T ITC-FR  FCOUI' T,rQnTIMT(2) ,CLI (2) ,POT(2) 

'■•iTrorP  FOATRI'Y,  -Oi.'‘FL,FT»JISH,'^PNTPAP,PIXYFUL-(2) 

TUTP-.FO  FOlXY,pP'. '’TF5:,FPi.!^MSH,''PLLN<!M,FPUHRS1,FPLnASN,FPLAPO 

pFAL  OT.lNTIMc  (1, 

XYT,IXY,< Y.PFXY'.OPXY 
APTtAF.A’O.oncs.praoo 
OOMMON  PFPOFS 

-OKMOM  ln,un'fu,ypf,no,unc 

common  FCONInT,  p D4TTM'',rL0,  FLTt  ^OT 

FOATPPY,  PPUFFL.FINTSM,ff’P»4TP4P  .FIXYFUL 
TO-Im-ih  FPLXY.FP.  “OPS  ,FPLP‘*^*',p“LLN3M,FPLH9SM|  FPlOASH,  FPUAPO 
OOMM-U!  rT.LKTIXF , OA YT,TYMP 
'•nM^nN  PLOTxtP(t  TO,  rJO) 

TF(FpLXY*FPLAPO+FPlPOFS,<;p.YFS)  30  TO  510 
rF(PPLP>IS1.EO.YFS)  ro  TO  511 
TP(F''LLNS1.F0.YP5)  r,n  to  512 
IF<PnLHR5>l.EO.YF5 ) 00  '^O  51T 
TP(F"L0ASH.Fn.YF3)  30  To  51<, 

510  ’'■(FPLYY.EO.YFS)  CALL  OTSXY 

. I''(F’'LAPO.EO.Ye3)  CALL  OISAOO 
TFfFPLPOFS.EO. YF5)  CALL  OISPOPS 
pCTiiPN 

511  CILL  OTSFMSI 
pPTU’N 

512  CALL  OTSLNSI 
PPTIJ^H 

51T  call  PfSHRSM 

PPT|J-',| 

511,  "ALL  OTS^ASM 

prTi,5„ 

* r»|n 


r i 
r' 
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o o 


■?unpoi)TiNE  oi5xy 


OrSPL4r  OP  SIGNAL  SP4CP  OATA 


j 

cn**aLPx  xvT(iP»,rxv(i2a,?»  ,xy(i’o,2)  ,PPxyT{i6>  ,ppxy(i?o) 

OOMPLPX  APT(lt),lP(12a,?)  ,APn(  130,3)  .POPS  (120, 2)  ,PPAPri(120) 
OOXPI  PX  PPPOPE (I’O) 
intpofo  ln,lnsel, yes ,mo,ln'; 

T-irc-.P7  FCO^^r^T,  P0ATr’‘(T(2)  ,PL1(’I  ,PLI  (2)  ,P0T(2) 

. ‘'■"■OFP  FOATRoy,  P0LSEL,PI*IISH,pfntpAP  ,FlxyFUL  (2) 

T'ITPOFP  FOLXy  ,P0'.  o 1PS,FPLP>1'm,pdllnSM,  FPL‘(‘?SK,FPLnASM,FPLAPO 
ppAL  nT,L>(TIMF  (2) 

30“HriH  XYT,lyy,<Y,FPXy’’,PF»Y 

COMMOU  APT,Ar,A3n,PpPS,PPa“o 

“OMMOM  pfpOFE 

:ommon  LN,LH''EL,ypp,No,LNr: 

:0«Mnn  FC0NIf'T,P05TTNT,PL0,Pl,T,PnT 

OO^HON  FOATP[;Y,pplSEL,FTNISp,pFSTPAR,PIXYPUL 

COOH-'AJ  FPLXY,  FP.  DOFS  , FPL P>*S>',P“L LM3M  .FPLHPSH,  PPLO ASM, FPL APO 
-OMMON  DT.LMIM-  ,DaYT,TY«P 

PPAL  MaPXY(12'i,l’9) 

'OMMOW  MADXY 

OATA  TARCFT/IHT/,  STAP/lH»/,Mr'0/H-/,  VEP/1HI/,9LANK/1H  / 

SFT  MAoyy  plank 

00  It  r=l,129 

00  10  J=l,129 

MAPyy (r, J)=9LAN< 

13 

'.ON'TNUE 

11 

TOnTthije 

C 

Spy  MAoyy  hop  GRID 

•00  13  1=1,129 

00  10  J=l,129,33 

HAPXV  (I, J)=HOR 

t?  OOM'TNUf 


c "‘^T  naovy  Yt”  r.pro 

TO  1'  T^i.irq.ip 
00  14  J=l,129 
mAOvy  (I,  J)=(/FP 

14  "ONTTMIIE 

15  •lON-'TNIie 

c STT  'ABOYY  VER  RCALE 
. 00  1’  1=1,129,32 
00  15  J=1,120,A 
M50XY  (I, J)=HO? 
15  I'OArTNUF 

I’ 

r,  SP'  xtOYY  HOP  scale 

00  11  T=l,l3y,j 

00  1 A J= 1 , 120 , 33 

•‘'PXVd,  J)rVP^ 


\, 


*' 


111 


ts  cofriNUP 

19  OONT'^HIJE 

C '’AT*  STAR*;  TO  *1APXY 

'^0  T = 1,1?0 
orxv (T)  = (PFXY  (I)  ) /2 
TXs'‘FAl.(PCXY(I)» 
TYsSTf'AGCOFXYdl  ) 

■^APvy  (i5»TX,65-IY)  =STAR 

20  ';nN’’TM',)S 

0 A-S'^'fl  "A^GFTS  to  HAPXY 
">0  ’I  T = l,  16  ■ 

•’'■XYT  (I  I r (PFXYT(  T»  > /2 
TX  = 'FAl (OFXYT  (I)  ) 

IV  = .'.T.iftr,(PCXYT(:  ) ) 

MA-’YY  (h6*IX,6;-IY)  xTAPF.ET 


21 

notr  IMUE 

r. 

rjFOT  Tn  TOP  OF  NEXT  PAGE  ARO  PPIHT  TITLE 

; shift 

TO  8 WEP  LN/INCH 

PPIR*  90  ,LNSEL ,D» YT 

, TY1F 

50 

PQPH  T (••l","CUR’- NT 
1 TIi-|",A10) 

oata  signal  space; 

LINEI", 

12, “ ; 9ATEI~,A10,“; 

oojM-  T5 

55 

FOPli-  (••!••> 

OPIN''  'O 

20 

PPPIA  T ( 3X//2X) 

„ 

C 

PPINT  GpA=H  OISXY  tip 

SCALE 

PPTN’90 

90 

<•  'P1ATI2X,”  -128 

-112  -96 

-90 

-64  -48  -32 

1 -if>  0 «^16  4’T2  *fth  *6^  »96 

2 ♦112  ♦128") 

C OPIMT  r,-ACH 
*1A?i<=7 


-in  ICO  HLN=1,129 

I'-(iAf'x.E0.7) 

so 

TO 

59 

IPC  'LN.SO.hFI 

Gl 

TO 

96 

TP  l'<LN.F0.i,6) 
G6  -n  9'J 

G") 

TO 

5’ 

0 ''I'IMT  wtth  vFF,  sriL-  nuinrps 
99  MVSr  1 

PPI  IT  61,HVS,  (HAPXY (I,MLN» ,T  = l,129) 

61  fOO  •AT(3X,I't*l21Al) 

•<A  = / = 0 

sn  ■’0  100  -• 

r O!  .NT  MITM  v'^R  AXIS  LAIFL  "Y" 

99  '•PINT  62,  (mAPXYI  I.hlNI  ,r  = l,l->1) 

62  Fno'iA''-(2X,lHY,TX  ,1TSA1) 

«apk=map<ti 
SO  *0  100 

2 "PTMT  HTTP  YFF  AXIS  LAO*-!  "AxxS" 

-7  PPT'I’  ft-*,  (NA'’XY(  1,11  H)  ,1  = 1, 1F2) 


RrVi  AxiiiiAHiF  (i'lK/ 
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ftS  F0l>M«T(2X,4HftXI5,lX,129Al) 

MASKsIARX*! 

r.O  ’O  luO 

C nqimr  PLAIN  MAPXY 

S9  PRI'iT  64,  (MAPXYd.MLN)  ,T  = 1,1?9) 

64  POP 'AT (7X, 129A1) 

MAPXrMARKTl 

r.O  *0  100 

100 

r;  e^I«lT  SCALE  AND  AXIS  LASPL 

=?TN'  60 
p^itr  65 

65  c’ORI.'M  ("S") 

’IIM'  S6 

66  nOMlTC*  “,10X,“X  AXIS~) 

C CLFA'’  ni^xy  FLAG  ”FPLXr“ 

FOLXV=N0 

RFTtPN 

run 


} 

tel  AVAllASIE  COP'/ 
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«;(lflPOUTINE  0XS4P3 


c ni'JPLAY  ''T  amplitude  and  phase  TEVTATION 

fOMPL^X  XYTdf.)  ,IXY<  130,2)  ,»Y  (120 ,2)  ,PFXYT  (16)  ,PFXY(120) 
COMPLEX  APT (If. ), A 0(120, 2) , APT (120 ,2)  ,onFS (13  0,2) ,PFAPO  (120) 

'OHOi.FX  PFoopt;  (i->0) 

X»(T‘^'T0  LN.LNSFL,  YFF,Nn,LM(' 

X.|Tr-Fi?  FCOMIHT,- OftTI'IT  (2 ) , r|  Of  2 ) , FL I ( 2)  , F TT  (2 ) 

’ 'TrOFP  FnATP'Y,  FPLSFL,FINISM,rr,|TP40,FTXYF|)U(2) 

r.iTO'TO  poLXY.FO.  P3FS,FPt.FHSM,FPLLN3M,FPLH<?SH,FPLOASM,FOLAPO 

o''AL  OT,LH''lMr  (2) 

XYT , I X Y, < Y ,oFXVT, apy Y 

COMMOM  40T,AC ,A>n,POFS,OFaon 
on  1MOH  pppop^ 

OOHMOM  LH.LN' e.  , Y"S,tin,LHC 

OO  frtOH  FCONIf'T,- n4’-I*IT,FLO,  FLT,FnT 

OOMMOH  FOATPOY,  0 0LSFL,FI»)TSH,rF4TP4P,FIXYFIJL 

■OMHOM  FPLXY,FP.oop5,poloHSH,F3lLNSH,FPLHPSH,FPLOASM,FPLAPO 

"OHHIh  DT,LN*I>i:  ,.nAYT,TYHr 

ofjl  MAPA0n(l?2, 121)  • 

'"OHMOH  MAPAOO 

OA'^A  AMP/ihA/,oh\ /lHP/,H02/lH-/,VF9/lH?/,RLAf4K/lH  /,LO/2HLO/ 
OAT.*  71,/ihL/,  ti/ihI/,  Tf)/HN/,o/l-(E/,0/lHn/,ll/lHU/,T/lHT/ 

C SFT  MAP.Aon  OlANX 

no  1 1=1,122 
no  •>  j=i,i2i 

HAOAoOd,  J)=')LAN< 

2 . nnn'-TnuE 

1 •n^nrnuE 

on  2 T=l,122 
nn  . 1=1,121,12 

C SFT  HAPa-'O  MOO  ORIO 

HAPA'’0(r,j)=H0? 

4 rOH'TflUF 

3 OnV'nirf 

C SFT  MAPAon  VFR  GRIO 
nn  F J=i,i2i 
on  1=1,61,10 
HAPAPOd,  J)=VER 

6 nnu’-j.iup 

no  ’ T=62,12?,13 
MAPAOOd,  J)  =VF^ 
nOMTTIIUF 

C ASSIGM  Apn  nATA  TO  MAPAPn 
no  A T=l,120 

TF( ’FALfPFAPOd)  ) .GF.28.)  GO  TO  10 
Ar^TAL  (OTAPOd)  ) 
r=' ’ tAG(PCAon(I ) ) 

TA?' 
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IPs" 

•1AP!VPn(IA»31,I)  = A*1P 
MAPAPOJIP+gZfDsPHA 
^0  '0  8 

10  MAPf.''n(58f  I)=7L 
•iAO'iPo<«;q,i)=7r 
MAP‘'P''t60,I)=ZN 
*<AniPn(51,I)- E 
lAP' 00(62, I)=nLlNK 
‘iRPiOO(63,I)sO 
riAPAO0(64,I)=n 
MAPAP0(o6,I)=T 

8 ':;onT’'N'iE 

C ejl^OT  TO  TOP  OP  NEXT  PAGE  ANO  P?I‘<T  TITLE?  SHIFT  TO  8 WER  UN/INCH 
P’lNT  P.O,LHSEL,01VT,TYMF 

SO  FO  = HA’' (“f.-ClIR^-MT  DATA  AHPLTTUOE  ANO  PHASE  OEVIATIONS?  LlNEI“,r2 
1,"  ? nATEI",A10,'*:  TIMFf.AlC) 
ppiRr  55 

55  P0OHAT(“T") 

C PPINT  GOAPH  legend 
PPI»r  GO 

60  COOHA- ("O",”  AsAHPLITUOE  OPtflATION?  P»phASE  OEVIATIOM”) 

C PRINT  TO'’  'CALING  INFORMftTION 


““IN’’  65 

65 

F3PaiT(“0”,10X, “SCALED  WITH  HOOPN  X.T  WALOES**) 

nc  T N-  T 0 

70 

FORNAT(“  “,qx,"-70  -20  -10  0 

♦ 10 

123  ♦33*1-70  -20  -13  0 OEG.  *10  *20 

e ♦•'C”) 

C PRINT  OPAPH  no  LOOP 
MAPXrll 
TVV=-f 
■>0  3 

irC'A^K.EO.ll?  ",0  TO  101 
!»■(  'LH.F0.50)  33  TO  102 
IFC'L’I.FO.Sl)  53  TO  107 

C PPI'IT  "LAIN  GRAPH 

"“TNT  7r, (NAPA03(I,MLN),T*1,1?2) 

75  FOP  iA’(llX,122Al> 

■<JP<  = NAP<^1 
00  '0  <3 

0 PPiM*-  -.P.IPH  WITH  SIDE  SCALP 
131  (ri/:lV(/*5 

“"TNT  3D,ITV,  (NAPAPOd.NLN)  ,T=  1,122) 

63  FOP  •AT(nx,I2,12.?A1) 

MJPVsO 
0"  -0  R 


miNT  rr.iPH  hiTm  Sin-"  AXIS  LA1PL 


O 


102  opit^T  OS,  (MAPAP1  (t,  HLN)  ,Ist,l?2l 

85  F0BMAT(1X,"TIME  IN" , 3X , 122A1 » 

NARX=NARK»1 
r.O  '0  9 

103  P”TNT  90f (NAPAPO d.MLN) , t=l ,122) 

90  roP'IAT(lX,"mLLtSEC.",lX,12?Al) 

NAOi^rHAfiK*! 

50  '0  9 

fn"  point  loop 

jNTTnue 

C oropiNT  TOTTCH  SCALING  INFOPMATTON 

P2INT  •'0 
’’IN-^  55 

c rffoin-'  lfgeno  at  bottom 
»RINT  fiO 

C PFTN2N  TO  STANOARO  LINE  SPAFING 
“’IN’  05 

95  F00MaTC"S“» 

c clfao  oi-olay  apo  flag 

FPL Aon=NO 

’■T!|0N 

r .'O 


Bbi.Av/iiLABLE  CUr 
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SnUPnuTIKF  OISPD^S 

c 

C OISPLOy  OF  PHASE  OEVIAnON  F!?EOllE-ry  SPECTRUM 

CO-PI.EX  XYT(16)  ,TXY(  l?a,?)  f XYd’O  ,2>  ,PFXYT  (16)  ,PFXY(1?0) 

COMPLEX  APTLIS), 1 P<1?0,?) ,AoO(120*2) f POFS (120,2) ,PFAPO (120) 
complex  PFPCF':(1?0) 

IUTF-pr  LN.LNSCL, YES,Mn,LMC 

TUT-'-.Fp  FC'1NI?.T,-0ATTNT(2)  ,e’LO(2)  , CL  I ( 2)  , ^OT  (2 ) 

’ TF-.FR  FnATK3Y,cP;.c-cL,rijr5M,FEYTPAR,FrxVFUL(2) 

FPLXY.FO.o^fs.fplc-Sm.FPlLUSM  ,FP.HRSH,FPLf)ASM,FPLAPO 
ocAL  M', LMTIMF (») 

CO“MCH  XYT,IXy,XY,PFXYT,PCXY 
CTiMCw  APT,ap,i3o,por$,pc4pn 
COIMOM  PFPOFF 
COmmim  LN.LN'^fL,  YCSfUn.LNC 
common  FC ONI NT, cost  I NT, FLO, FL T , cot 

common  foathov, colsfl,fintcm,cfmtoar,fi<yful 

common  FOLXr,FP^PJFC,FOLCMSM,Fpt.LNSM,FPLHRSM,FPLOASM,FPLAPO 

COMMON  nT,LFTIM",OAYT,TYMF 

-OMMON  PLOTMflP( ITC, 133) 

OATA  STAF,/1H*/,YFR/1HI/,H0P/1M-/,PLANK/1H  / 

C SET  OLOTMAP  RLAMK 
00  1 T=l,l?l 
on  2 1=1,121 

PLOTMAP(I,J)=3. ANY 

2 CON'IHUf 

1 CCNTTNUE 

C SET  HO’  OLCTMAP  GRID 
•00  3 1=1,121 

on  - 1=1,121,12 

PLOTMAPd,  J)sH3R 

4 CONTTHUE 

3 CONTTN'IF 

C SET  VC"'  PLCTMAP  GRIO 

00  = t=i,121,?o 

OO  <■<  J=l,121 
PLOTHAP(I,J)=YiR 

6 continue 

5 CONTTUUE 

C SET  POcr  plot  DATA  INTO  plOTmAP 
no  7 J=l,64 
CMAO=CARS (PFPOFJ ( J) ) 

•'AO:'-Msr/{500./l20.) 

:F(  16C.r,F.12C)  MAGsllO 
YsO.OTTC* 1 

PLO^MSPIX*?,!  21.M4';)  sSTAR 

7 CONTINUE 

C FJEC’-  -n  TOP  OF  NF<T  PAGE  ANO  P'^INT  TITLE!  SHIFT  TO  « WER  LN/INCH 

c.D^U*  na,LNGFL,P\VT,Ty'(r 


AVAIIASIE  COi' 
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50  FOR*1»T(“l-,“CUR^"NT  OAT#  ohAFc  OEVIATION  FREOUPNCV  SPECTROH;  LINE 

; OATEt“,AlJ,“t  TIMFI“,A10) 

PRINT  <55 

55  FOOHAH“T“> 

C PRINT  TOP  LAPEL 
PPINT  60 

60  FORHnT(“0”»15X,“rRFOUENCy  IN  MFRT7“) 

C PRINT  TO“  scale 
' INT  65 

65  (~0",15X,"a“,16X,"200",17X,"400“,irx,“600",17X,“800",16X,“10 

100",15X,“12J0") 

C PRINT  '.PAPH  LOOP,  SIDE  SCALE,  ANO  LAREL 
MAOXrll 
ITV=-50 

00  fl  MLN=1,1?1 
IP(HARK.EO.ll)  30  TO  21 
ir(MLN.Fa.54>  50  TO  24 
ir(HLN.E0.53)  GO  TO  23 
IFC'ln.F0.52>  50  TO  22 

c pptnt  "lain  graph 

“PINT  7i.,(PL0THlP(I,HLN)  ,I  = l,12l> 

75  F0RHA’«16X,121Al) 

HAP<=HAP<,1 
GO  TO  8 

C PRINT  r.ojPH  WITH  3I0E  LABEL 

22  tPTNT  80, (OLOTHAS (I, MLN) ,1=1,1211 

80  FORH'.T  (IX, “FREQUENCY”, 6X,  121A1) 

HAPTsHARK*! 

. GO  ’0  8 


23  ORTN’  .35,  (PL0TMAR(I,MLN)iI  = l,l’H 

85  EOphAT  (1X,"SFE’:TPU*'",TX,  121  a 1) 

HAOi'rHAPK,! 

GO  ’0  8 

24  “PINT  90, (PL0TH1P(I,MLN) ,1=1,121) 

90  PQP'^A’  (lX,”MAGNt  TUnE“,6X,121Al) 

HA“<=HARK»1 
GO  TO  8 

C PRINT  G“APH  WITH  SIDE  AXIS  SCAL'^ 

21  IYV=IYV-5a 

• POINT  q5,IYW, (PL0TH#P(I,MLN) ,T=1,1211 
95  ■•noNA*(iX,llX,I,,121Al) 

lAPt'sO 
GO  TO  8 

C “NO  GRAPH  LOOP 

8 CONTTfiuE 

C npoor.iT  OOTTON  scale 
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“RTIT  65 


C BFPPIN’’  BOTTOH  LAB'L 
Pi^INT  60 

C PFTIIPN  TO  BTANOARO  LINF  SPACI*<G 
PRINT  100 

too  FO?N('T(~S~» 

C CLFAP  OTSPLAY  FLAG 
r^LPOFS=NO 

P^TU’H 

FNO 


5UHO0IITIME  0ISW5H 


“0>1PLFX  XVT<lf,),rXT«12l),2)  ,«T(120,2»  ,PFXFT(16)  ,PFXY(12  0) 

nOMPLFX  *‘>T(16),4P(120,2>  ,«»0  (120,21  ,P0FS(  120,2)  ,PFAP0  (I  20» 

OOMPt.FX  FFPriFS  (1?0) 

IMTFrcp  lN,LH3FU,  V'-S,N'),L‘('' 

lUTc-.r's  FC0NI‘(T,Fr)ATI‘|T(2)  ,flO(2)  ,FUI  (2)  ,f0T(2) 

imr-rn  niTFOY,  FPLS<-L,FIMT*M,=’rsTP4P,FIi<yFijL  (2) 

• ^rrro  FPLYV,FP.o-IFS,cpLFM<;M,rpLLNS*1,FPLHPS1,FPtO»SM,FPL»PO 
■'4L  DT.LNTI’if  (2l 
OOMH-'M  XYT,:yy,<Y,ppxYT,PFXY 

^O-oiTN  APT, -3, pdfs, PF#*" 

OFooc*; 

'rir^MTN  LN,LN‘a,Y':S,nr),L‘)'' 

"OH.*  1,1  FCn‘4l>iT,Tf)4TTNT,rLn,FLt  .TOT 

:0-<«'’N  FTlTPr!Y,FPL':FL,FINIFH,FC,,TPAR,FlXVFlJL 

FPLXy,FP.oiF«;,FPLFMS“,FOLLNSM,FPLHRSH,FPLOASM,FPtAPO 

0T,LHT  liF  , jAYT.ry^lc 
:o'>nN  OLOT-aPdxj.isj) 

■>FT||->,( 

r.jo 


120 


c 


SlinPOIITINE  0ISLN5H 


iomplfx  yyT(i6),rxY{ 120,2) ,xy(i20,2) ,ppxyt(16) ,pfxv(120) 

OOMPUFX  APT(lf>),4P(120,2)  .APT  (120, 2)  , PDFS  (120, 2)  ,PFAPO(120» 
OOMPl.FX  PFO0FS(l?0) 

LH,LNSEL,  Y‘='';  ,IIO,LUf: 

TNT'^'-F?  FC0NIMT,p04TT>|T(2)  ,rLO(2)  ,FLI  (2)  ,FTT  (2) 

TMTFnFD  FnATP()Y,FPL':rL,FINT<?M,FENTPAR,FIXYFUL(2) 

T ■TF''.F?  FPLXY,FP.PlFS,PPLFM5M,FPLLrl3*(,FPL‘(PS’l,FPLOASM,FPLAPn 

’"5L  n’.LUTlHr (2) 

XYT,lyY,<Y,PryYT,PFXY 

“OHHPM  APTtAF.A’OjPpPS.PPAPn 

pfoofs 

OTH-nn  LU.LMSFL  , Y"^*"  ,N0,LN; 

COHPPN  FCnniNT,  =’OJTIMT,'"LO,  FLYt^OT 

OnM’IPN  FOATR^y,  '•PL':'’L,FIMTFH,'‘P'(TPAP,FIXYFUL 

OOHHIH  FPLXY,FP.PlF5,rFLFtY3M,PPLLNSM,FPLHRSM,  FPLDASH.FPLAPO 
CO'IMON  OT,UNTIH-:  ,0A  YT,’Y>*F 
COH1PN  PLOTMAPdSO,  130) 

PFTiPH 


r. 


F'iRpnuTiNE  OISOA3M 


-OMOLrx  XYT (Ifc), I XY( 120,2), XY (120,2) ,PFXYT( 16) ,PFXY(120) 

COMP'  PX  APT(16) ,1”( 120,2) , A 00(1 2C,2) , POFS ( 120 , 2) , pfAPO (1 20) 
“OMOlfx  PFP0FS(120) 
iNTFopo  ln,ln<;£l, ycs,No,LM: 

t'itf-.Po  FCOSIOT,"OSTTf|T(2)  ,rLl(2)  , pl T ( 2)  , p OT  ( 2 ) 

lOTP-.ro  FOATRnY,FPLPFL,FtflI3M,FP«jTPAi?,PIXYFUL(2) 

! Tr-r;?  FPLXY,FS.O0FS,>'OLF“3“,polLN5M,FPL4PSM,FPL0ASM,FPLAP0 
■AL  OT.t.NTI.'IF  (2) 

'OMM-Ifl  XYT,  lYY,  < Y,ocyyT,  PPXY 
AOT,AP,A->0,opFT,PP4Pn 
' pcpir-i 

'OMHOU  L0,LH''FL,  V£5,NO,LM" 

: T'MT»(  f coMi'iT,  p.o'iTr*(r,FLn,Fi.r  ,por 

•;  OM'dO  F047  P"'y,FOLPPL  ,FI')I5-t,FPSTOAR,FIXYFlJL 

-0‘<M''H  FPU»Y,F0.  o ->£■?,  pflP'"7m,  FPL  LN3M,  FPL  HPSM,FPLO  ASM,  FPL  API) 

''V‘’nn  nT,LN’'IM',0AYT,TYxP 

'0  '«-'N  PLOT'110(  1 Y j,  I 30) 
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rso  NOS/BE  Lt.l4T  ryBP  CMPS  08/31/77 

14. -JE  . FPOM  /fa 

14.  yf  .'•  ^.la  ooOOBrf,^  W053S  - FTlF  IMPUT  , OC  10 
14.  ■»<  .'^’.MIs,r'^lJ1Cim,''T''f1.  r-'7  0?6Bt'^T^'T1*J''P»''«c-7 

14.  I-'. . ' ./T,  u: 

14. 7f,  .CTN  ,'’  = P'IMrHB  .flPTi  '.Pr  0 . 
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